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ABSTRACT 

In this study a linear two-level baroclinic model is 
formulated in which the zonal steady-state current is 
perturbed by surface friction and surface sensible heat 
exchange. The perturbation quantities are assumed to be 
either independent of latitude or varying sinusoidally with 
latitude. 

Investigation of the phase velocity response of the 
model indicates that sensible heat exchange decreases the 
amount of instability of short and intermediate waves, but 
increases the longwave instability. The combined influence 
of friction and diabatic heating is to create unstable iong 
waves which travel eastwards at reduced speeds and which may 
be stationary or retrogressive in nature depending on the 
speed of the mean zonal current. 

Phase and amplitude characteristics of the existing 
perturbation waves suggest that the thermal wave necessarily 
lags the geopotential wave during amplification while the 
static stability wave prefers a lead position. In addition, 
maximum downward vertical velocity and maximum surface 
sensible heat transfer to the atmosphere occur between the 
ridge and downstream trough position. 

It has been suggested that blocking ridges may be a 
form of longwave baroclinic instability and that sensible 
heat transfer plays an important role in their formation. 


Comparison of numerical predictions of unstable longwave 
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amplificaticn with observed North Pacific blocking space and 
time scales indicates that without the presence of sensible 
heating in the model the development of blocking waveforms 
is not feasible. With sensible heat exchange stationary 
unstable waveforms are obtained which have growth times and 
wavelengths that are similar to observed North Pacific 


blocks. 
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CHAPTER 1 
INTRODUCTION 


1.1 Theories on Blocking-Ridge Formation 

Blocking-ridge activity is one of the most visible 
forms of longwave growth in the atmosphere and has long been 
recognized as being a major influence on regional climate. 
Such enduring ridge formations, which are found to be more 
frequent and pronounced during the winter, can displace 
passing synoptic disturbances considerabie distances to the 
north or south of normal storm tracks. This then invariably 
leads to anomalies in such parameters as surface temperature 
and precipitation when compared to the non-blocking regime. 
Although the climatic effects and behavior of blocking 
ridges have been well documented, the dynamic causes 
surrounding their initial formation and development are 
still not fully understood. 

In an early consideration of blocking phenomena, Yeh 
(1949) suggested that the transfer of energy by dispersive 
waves in the atmosphere might explain the occurrence of 
sustained ridging. These waves are created by a process of 
rapid adjustment between the pressure and velocity fields 
and energy propagation occurs at a group velocity which can 


be very different from the speed of the prevailing westerly 
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current. By examining the mia hover tereee of a solitary 
wave within an incompressible uniform atmosphere, it was 
Shown that cnce a pressure rise or fall was created near the 
pole it would remain there without being dispersed. At 
latitudes other than the pole the solitary pressure wave was 
found to retrogress while slowly spreading out and 
dispersing, a characteristic often exhibited by atmospheric 
blocks. From these results it was concluded that blocking 
activity was likely a high-latitude feature of atmospheric 
flow which has greater intensity, a longer lifetime and 
moves with slower velocities as latitude increases. 

Although this concept can describe a few of the 
prominent properties of the blocking process it represents 
others inaccurately. For instance it overpredicts both the 
extent of the motion and the rate of dispersion of blocking 
ridges at lower and middle latitudes. Figure 1.1 
illustrates the behavior of a solitary dispersive wave at 
4ON and 7ON. The curve representing the theoretical results 
for 40ON shows that after a relatively brief period of 32 
hours the wave amplitude has dropped to 1/3 of its initial 
value and has moved upstream more than 6000 km. At 70N 
after two days the same solitary wave had retrogressed at an 
average speed of 880 km/day and its amplitude had diminished 
by approximately 1/6 of its original value. The 
Characteristics of the wave at 40N are not representative of 
actual blocking behavior, while the wave speed and 


dispersion rate at 70N are more indicative of the last 
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Fig. 1.1. Dispersion of an initial solitary wave. The dashed 

curve is the initial waveform at time t = 0; curve (i) is the wave 
profile at latitude 40N after approximately 32 hours and curve (ii) 
is the wave profile as seen at 7ON after approximately 2 days. The 
abscissa is in units of 2200 km, while the ordinate is in arbitrary 
non-dimensional units of amplitude. The rate of dispersion of the 
wave amplitude at both latitudes increased significantly as time 
progressed. The speed of the retrogression at 70N remained constant, 
however, accelerated retrogression was observed at LON as time went 


on. Curves are taken from Yeh (1949). 
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stages of blocking when eiteidye is deteriorating. In 
general the theory does not describe satisfactorily the 
development or persistence aspects of the phenomenon nor 
does it explain why atmospheric blocks tend to show a 
longitudinal and seasonal preference in their formation. 
Other investigators such as Berggren, Bolin and Rossby 
(1949) have hypothesized that blocking in the form of a 
wave-like perturbation could be analogous to that of a 
hydraulic jump in open channel flow. Continuing this 
approach, Rex (1950a) employed a hydraulic analogue model as 
devised by Rossby (1950) and a momentum transfer argument to 


arrive at a critical velocity given by 


ae 
FE TESS iets 


where 4 is the linear Beta-plane approximation to the north- 
south variation in the Coriolis parameter f 

(i.e. fiy)=f,+ yi. and a is the half-width of the current 
flow. Whenever the actual zonal current speed exceeds u 

two permissible modes of flow exist with one being at a 
higher energy level and velocity than the other; a necessary 
condition for an analogy to a hydraulic jump. From (1.1.1) 
it is evident that decreasing the stream width and 
increasing its flow speed meant lowering the critical 
threshold of the flow such that an external impulse to the 
stream might cause the jump to a higher energy state. 
Therefore a region aloft with a narrow strong jet and 


anomalously low thickness values (shallow stream depth) 
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would lend itself more easily to the formation of a blocking 
ridge. The propagation of the atmospheric ‘jump' once it 


has formed was determined to be 


c= V-n(n+1) (1. 1. 2) 


where N=u/u, is the critical velocity ratio, u is the 
immediate upstream current velocity, and n=a, fa is the 
expansion ratio of the upstream and downstream current half- 
widths as indicated in Figure 1.2. 

After applying the theory to a number of test cases, 
Rex (1950a) concluded that no major contradictions were 
apparent and that the concept gave credible results in 
predicting when blocking might occur. However, one argument 
made against the hydraulic jump theory is that it lacks a 
feasible energy source required to excite and then maintain 
the developed wave state associated with occurrences of 
blocking in the atmosphere. There is also some question as 
to how applicable the idea of a finite current width is to 
the real atmosphere. In narrow ‘well-defined stream widths 
when the term a is small, calculations of u are apt to be 
unrealistically small and the resulting inordinately large 
values of N and n do not lead to the correct interpretation 
of blocking behavior. In situations when the flow 
boundaries are not so well-defined there is often some 
arbitrariness in the selection of the half-width and 
therefore some uncertainty about the value of the 


predictions. 
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Fig. 1.2. The simple two-dimensional steady current envisioned in 
the hydraulic jump analogue model used by Rex (1950a). Zonal current 
’ speed u is assumed constant over the cross-section A, as is u, over 
cross-section B. The ‘'jump'' can occur according to the theory if 

N = u/u. > 1, where the critical velocity is defined as ae Ba*/3. 
Current flow is assumed to vanish (U = 0) in the environment 


surrounding the channel. 


In an investigation of the zonal index cycle, Namias 
(1950) alluded to the possiblity that blocking action was a 
manifestation of baroclinic instability in the atmosphere. 


Later, after studying the interactions of an ocean- 


atmosphere environment, Namias (1959) proposed that the 


thermal energy anomalies contained in the upper ocean may 
provide the necessary energy source for occurrences of mid- 


ocean blocking. Other studies on searatmosphere heat 


rn | earn. 
\ he 1 aa 
re iy if 
ech ri 
) wi ee ae eS edie eat | ei § ae Wi ie 
st { / 
Size! Beit as st 9 Magi awe S008 ed ost ” a ai Co i i a es rn | 46 AS 
. 3 2 © : a). - c ) f £, 
= 
“* J. 7 a + ng te 
i 
: r oo . + 
tae oye een) EEE RTs ia te iG a “Re 
‘ . =) owt : 
- « = Wes he bac 4 , ' a 7 acs i 
= Satee> eagle Geen 7) 
a) er p 1 
t veh i om 4b 
5 i ky } 
PU ; a Bs Nt: 
; jh i en 
= 
. y = ; 
se ah dee ad — ee! = me Ge 
. ‘ie 
> a 4 calle J iM 
Sewn clone ie yt 
~ 
\ 
daha hh al vai To a shite 
F 
‘ i Fi} “ih "4 f\.9 mr 
= Le Rye i 
. , . i ‘eAbes Aa ie © bh ‘ 
eT v3 ? A tO) Do ieee ere be wid 


hee | ae ‘ - 8A 0 
user P bath Hy attr noaed ‘noe Vea | 


a re 2 


sore ry: ‘ 7 iTeos 
7 a ee bond al ps (> 'analin Anata 
jase ian ey Sit ne tie ue iain oF Lee tnd 


- ir) ak ane r ‘i ae ie he 


nah te Z a : = it / ; 4 eh . . 
: ; :) ae _ A . ns “4 iad ola 


‘ An 


~ 


Ja 
rere ipsa’ meh ick ee pene econ al 
ve 

Rae An ksss par Acerde. 24e3 yiderween en Se RbOw kes PM a 


a) 


seipaqtiat # oe ete hate sads) “Fark: one, tw 7 
‘mower he 6. e088 rubas tas was Satatbeany abrtp eto 
MAAS. sng2, baconwodgh Pus): Fitank cpap eovaimnetn: | 
is aps hjique ana’ yf, papas 409" ae anon UpRene- ceaamay, ae 
big ehuaer3b: iti GO? "nase todayy Waksvavisna dir witmege’’ a 
a ae Kant Miner diespel malate a sete 7, 
7 hee 7 | : | 


i 
a 


transfers, however, have iaiehech that fluctuations in the 
ocean's surface temperature do not normally play an active 
role in governing thermal energy exchanges ona time scale 
such as is associated with the development of blocking 
ridges (Kraus and Morrison,1966; Clark,1967). It would seen 
that sensible heat transfer in response to surface air 
temperature fluctuations could form a more viable energy 
source having a suitable time scale of several days to 
several weeks. 

Theoretical results which have been interpreted as 
supporting the baroclinic instability concept of blocking- 
ridge formation came from work done by Haltiner (1967). 
Using a linearized two-level baroclinic model he has 
indicated that with the input of sensible heating from the 
surface a region of possible stationary or retrograde 
unstable wave growth is created for longer wavelengths in 
the range of 7000-9000 km. The instability of long waves in 
general was found to increase gradually as the amount of 
sensible heating was increased while waves on a near- 
synoptic scale became slightly more stable. These results 
were obtained for representative 50 kilopascal level winter 
values at mid-latitudes. 

Encouraged by this theoretical evidence, White and 
Clark (1975) examined monthly mean atmospheric data taken 
over the North Pacific covering a period of 1950+1970. 
Their study found that blocking action was most often 


observed over the North Pacific ocean at a longitude of 160W 
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to 1708 with a dominant space scale of 6000 to 7000 km. In 
addition they were to discover that the distribution of 
anomalous surface sensible heat flux was strongly somcinad 
with the anomalous 70 kPa heights associated with the 
development of blocking ridges. Below normal heating values 
were found to occur under the position of the ridging and 
above normal sensible heat exchange was associated with the 
adjacent trough regions on either side of the block. It was 
the conclusion of these investigators that sensible heat 
exchange between the ocean and the atmosphere was likely 
very important in the formation of blocking ridges and that 
the application of baroclinic instability theory to this 


type of phenomenon seemed plausible. 


1.2 Two-level Baroclinic Model 


wi 


Many of the baroclinic models which have been utilized 
in numerical prediction are of a simple two-level design. 
The procedure usually consists of applying a vorticity 
equation at two isobaric levels in the modelled atmosphere, 
(e.g. 25 kPa and 75 kPa) and a thermodynamic equation at an 
intermediate level, (50 kPa). The sum and the difference of 
the vorticity equations over the upper and lower information 
levels constitute averaged and thermal equations, 
respectively. The system of equations is then closed by an 
equation of state. Analysis of the dynamics of the model is 
possible by employing perturbation wave solutions with which 


the equations are reduced to that of the linear type and a 
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frequency equation is obtained. 

Early studies using this approach assumed that the 
motions in the atmosphere were both adiabatic sabe 
frictionless. Using a quasi-geostrophic linearized two- 
level mOdeis Holopainen (1961) was able to incorporate the 
effects of surface friction into his perturbation equations. 
For adiabatic frictionless flow, using representative 
atmospheric values, it had been established that the 
perturbation field contained a unstable component over a 
small range of relatively short wavelengths. The effect of 
introducing surface friction was to broaden this range 
slightly while reducing the component amplitude of the 
previously-existing unstable wavelengths. Haltiner and 
Calverly (1965) using a similar formulation investigated the 
influence of surface-based friction on the interactions 
between the averaged wave and the thermal wave. Their study 
gave evidence that the phase characteristics of the two 
waves were markedly affected by friction. Results showed 
that decreases in the limiting phase difference were 
predicted when the value of the surface drag coefficient was 
increased. 

Applying a different formulation of the hydrodynamic 
equations, Lorenz (1960) designed a two-level non- 
geostrophic model which allowed for full variability in the 
Coriolis parameter and the static stability. Gates (1961) 
has investigated the linearized form of this model under 


conditions of frictionless adiabatic flow. His results 
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suggested that the unstable perturbation component was 
significantly larger using a variable static stability than 
for cases in the usual quasi-geostrophic Poren iting when 
this stability parameter was held constant. The 
perturbation equations as derived in this study formed the 
basis for the previously-mentioned investigation into 
sensible heat influences on baroclinic instability by 
Haltiner (1967). 

The dynamics of baroclinic waves have not been 
investigated by analytical methods alone. Kikuchi 
(1969,1971) used spherical harmonics for the time 
integration of the vorticity and thermodynamic equations in 
a quasi-geostrophic two-level model. These numerical 
investigations were aimed directly at reproducing the 
occurrence of blocking action and at gaining a better 
understanding cf the effects of diabatic heating and 
mountain ranges on the behavior of such longwave growth. 
Included in the diabatic heating were the release of latent 
heat and radiational flux divergence as well as sensible 
heat exchange. Numerical computations showed that quasi- 
stationary longwave amplification similar to blocking 
processes occurred in the model whether orography was 
included or not, The main effect of the inclusion of 
mountains in the model was to lengthen slightly the duration 
of the computed ridging and to create preferred longitudinal 


regions where blocking occurred more frequently. 
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1.3 Intentions of This Study 

In this study a two-level baroclinic model which 
incorporates a variable static stability will be developed. 
Terms which will describe the input of surface sensible 
heating as well as lateral and vertical turbulent heat flux 
in the atmosphere are included. The effects of introducing 
surface friction at the lower boundary of the atmosphere and 
those friction terms arising from turbulent diffusion by 
small eddies within the current flow will also be examined. 

The governing equations will be linearized by 
neglecting squares and products of the perturbation 
variables. The dynamic features of the linear model will 
then be analyzed. The emphasis of this analysis will be on 
the stability of the baroclinic waves in response to varying 
values of sensible heating, friction and other variables as 
well as the interactions between perturbation waves. 

Of particular interest, in view of the findings of 
White and Clark (1975), will be the behavior of the so- 
called long wavelengths in the possibility of explaining 
some aspects of blocking-ridge development. It cannot be 
expected, however, that any perturbation method will be able 
to portray accurately all of the characteristics of the 
longwave amplification we observe in the real atmosphere 
such as those found in the blocking process. One feature 
that may be difficult to explain by perturbation analysis is 
the persistence of blocks. after having gone through a 


period of unstable growth, the ridge can be seen to remain 
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in a non-amplifying developed state for an extended length 
of time, implying that scme form of stable mode is 
accquired. It seems likely that this is beyond the 
prediction capability of a linearized model and, therefore, 
what will be of prime concern in this study is the period of 
time when instability is realized and wave amplitude 


increases to a developed state. 
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CHAPTER 2 


DEVELOPMENT OF THE MODEL 


In designing a simple baroclinic model some manner of 
modification is required to reduce the system of 
hydrodynamic equations to a less complicated form. In order 
to yield what would still be considered compatible results 
it is desirable that the reduced relations are consistent 
with any conservation principles expressed by the exact 
equations. The governing equations which describe 
atmospheric motions exhibit several such principles 
regarding the production of relative vorticity and energy on 
a global basis. They are: 


(a) the mean generation of relative vorticity 
over a global pressure surface is zero; 


(b) total energy under reversible adiabatic 
processes is conserved; and 


(c) the sum of the kinetic energy and the 
available potential energy is conserved. 


One method of simplifying a given equation is to 
compare the magnitudes of the various terms within it and 
then neglect those which are of smaller order. The orders 
of magnitude of individual terms are based on time and space 


scale considerations of the phenomenon to which the equation 
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is being applied. This approach is therefore called scale 
analysis and was introduced by Charney (1948). Application 
of the scale analysis technique while checking to maintain 
the integral vorticity and energy .ivariants allows the 
governing equations to be reduced in a consistent manner. 

A reduction in the complexity of the equations can also 
be made by using some form of general approximation 
pertaining to the atmospheric velocity field. The type of 
approximation used in many models is the quasi-geostrophic 
assumption in which the horizontal velocity field and 
vertical component of relative vorticity are evaluated by 
the non-divergent geostrophic wind. Charney (1947) has 
shown that assuming the horizontal velocity to be 
geostrophic has the effect of filtering out high-speed 
inertia-gravitational waves which occur in the equations of 
motion and thereby a significant simplification is achieved. 
The conventional quasi-geostrophic two-level baroclinic 
eodel (Charney and Phillips,1953), however, has inherent 
energetic inconsistencies which result in prediction errors. 
For example, there is a tendency for the model to over- 
predict the amount of kinetic energy in the atmosphere. 

Such inconsistencies and their effect on the model's 
performance are discussed by Gates (1961). 

A simplified set of dynamic equations which possesses 
certain consistent integral properties with respect to 
energy and vorticity has been formulated by Lorenz (1960). 


For a filtering approximation the quasi-geostrophic 
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assumption was replaced by a relation known as the balance 
equation. This approximation is derived from the following 


form of the divergence eguation 


Vd- Sareea Vie ee (2.11) 


where is the geopotential, V is the horizontal wind for 
which u is the west-east x-component and v the south-north 
Pecaupouenth 7 ai 7 )x + 5 Yyyed =0x¥ is the vertical 
component of relative vorticity, iy 4 and k are the unit 
vectors and J represents a Jacobian in the normal 
differential form. If a stream function Y ana a velocity 


potential X are defined such that 
a > ~ 
Yee Vy + Va (23122) 
V, « ri (2.1.3) 
i -4 VV (2.1. 3) 


—_> 
in which vy and vy are the rotational and divergent portions 


of the horizontal wind, then (2.1.1) takes the form 


Saad CAD ade 
Tye DPV) ea (ER) 2S oo ans 


Thus the balance equation as given by by (2.1.5) implies 


that a continual equilibrium exists between the rotational 
portion of the wind and the geopotential field. 

The prognostic equations used with (2.1.5) were the 
adiabatic thermodynamic equation and the vorticity equation 


neglecting friction. The former can be expressed as 
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where®=dp/dt represents vertical velocity, t is time and @ 

the potential temperature. Potential temperature at a given 

pressure p (kPa) and temperature T (°C) is given by the 


relation 
100 t Sccctes 
‘ @= ie / (2.1.7) 


in which K=R/cp is a dimensionless ratio of the specific gas 
constant of dry air R and the specific heat capacity of air 

at constant pressure Fai 
The version of the vorticity equation which was used 


can be written as 


ey. i Ulef)eayp (Lf) 7 V Vv x Fro O (2.1.8) 


By applying (2.1.2)-(2.1.4) relative vorticity can be 


related to the stream function by 


Se Rs Vy -V¥ (2.1.9) 


In addition if ¢d represents horizontal divergence VV, then 


Aon V. : an (2.1.10) 


Therefore, under hydrostatic conditions defined by 


ark 


where p is density, the continuity equation can be written 


4 
P (ate tt) 


in the following form 
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Employing (2.1.2) -(201.4), (2.129) and (2.1.12) the 


prognostic equations (2.1.6) and (2.1.8) can be written as 


2 T(¥0)+ 7697) +4 (v8) - 0 ae 


ROY TUT Vel) OX-O UL) sorb 9 
+ (Vbef)S + Yo VE 70 (2.1.14) 


Through a systematic analysis of the integral 
constraints, Lorenz (1960) has shown that together (2.1.5) 
and (2.1.14) have invariant energy and vorticity relations 
Similar to the exact divergence and vorticity equations. 
Thus, with the addition of (2.1.13), they form an 
energetically-consistent system of equations. 

This system of equations can be reduced further if the 
phenomenon to which the motions are being applied is of a 
evnoptic or larger scale. After (2.1.14) is subjected to 
such a scale analysis (cf Haltiner,1971) and smaller terms 
are omitted, a second approximation which maintains the 


integral vorticity constraint is 


S770 TY vi pit) + Vx-Vf- Ss aie (2. 1215) 


An equation originating from (2.1.5) which, when used with 


(2.1.15), satisfies the energy constraints is 


V¢-V#- VY-fv'}-0 (2.1.16) 
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Equation (2.1.16) is sometimes called the ‘linear’ balance 
equation. The set of equations (2.1.13), (2.1.15) and 
(2.1.16) will form the basis for the model equations to be 
used in the present investigation. 

It should be pointed out that for (b) and (c) to be 
true when hydrostatic conditions are assumed, tHe kinetic 
energy contained in the vertical component of the motion 
cannot be included in the overall total kinetic energy. 
However, it can be said that the amount neglected is an 
dnsignificant fraction of the total energy involved in 


synoptic and larger scale systems and, therefore, the 


omission is normally of small consequence. 


2.2 gensible Heating and Static Stability 
In order to incorporate sensible heat exchange in the 
modelling the First Law of thermodynamics will be used in 


the form 
do ; 7 
ol EF = H . (25 2511) 


where H is the rate at which external heat is added per unit 
Mass (i.e. H=dH/dt). By applying (2.1.7) the First Law 
states that the total rate of change cf potential 


temperature is 


ici 
a i F ide (262-2) 


We 
This is simply the diabatic form of (2.1.13). 
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The rate and direction of sensible heat flow across the 
interface of any two adjacent environments is largely 
dictated by the local difference in the respective 
temperatures on either side. In the atmospheric situation, 
where the air is in mction over a comparatively fixed 
underlying surface, the speed and turbulent state of the 
overriding air flow are also of primary import in 
determining a heat exchange rate. An estimate of the 
turbulent exchange of sensible heat at the earth's surface 


can be made using 


H. = 6 V, be Ze) (2.2.3) 


where T is the surface air temperature, T the underlying 
ground or sea temperature and V; the surface wind speed. 

The term C is determined from atmospheric parameters (such 
as the surface drag coefficient) and is constant for a given 
value of V.. Equation (2.2.3) is discussed in detail by 
Budyko (1956) and generally has been applied to calculations 
of turbulent heat transfers over water surfaces rather than 
over the continents. Haltiner (1967) has used this form of 
surface heat input in his model. As a further 
approximation, however, Haltiner (1967) assumed that the 
heating was distributed through the atmosphere according to 


a vertical profile of 
° -/p\r 
A(p) - AE) (2.2.4) 


for which H(£) represents the rate of sensible heating 
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occurring at pressure p, p, is surface pressure and r is a 
constant taken to equal 2 for his investigation. 

In the present study no assumptions will be made as to 
a predetermined vertical distribution for sensible heating. 
Rather, considerations of the turbulent nature of the 
atmosphere will be used to provide a mechanism for 
vertically transporting heat. The process of heat transport 
by turbulent eddies in the air is known as thermal turbulent 
diffusion. if H; denotes the sensible heating rate at some 


level j in the modelled atmosphere then it can be separated 


intc two parts 
Be Re br ae 


where R and S are the vertical and horizontal turbulent 
thermal diffusion terms, respectively. 
The vertical diffusion of sensible heat will be defined 


by 
dh, 
eas | 
K, i: (2. 2.6) 


in which g is the gravitational acceleration. Here hj 


represents the vertical heat flux given by 


hy 9% ky | (32) : r'| (e2n7) 


where K, is the vertical kinematic eddy diffusivity 
coefficient and I” (°c/kPa) acts aS a measure of the 
countergradient. The concept of an atmospheric 


countergradient flux is used to explain observations of the 
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occurrence of some upward heat flux even under slightly 
stable conditions. The lateral diffusion term is to be 


determined by 
2 F. x . 
Sask pee eesioll @ (202.8) 


where K, is the horizontal kinematic eddy diffusivity 
coefficient. 

To establish a two-level model which contains a self- 
determining variable static stability the thermodynamic and 
vorticity equations are to be applied at the same 
information levels, as opposed to the normal situation when 
the thermodynamic equation is used at one level only. 

Figure 2.1 shows the vertical structure of a two-level model 
where the isobaric surfaces 1, 2 and 3 represent the primary 
information levels. Therefore, with the thermodynamic 
equation employed at levels 1 and 3, a measure of the 
variable static stability at level 2 can be obtained by a 
eeEpical finite-difference representation of 967) p in the 


form 


(28) 2 8-2 , b-8 
fp Pa hy. Ps-Pz (2.2.9) 
From (2.2.7) it is clear that the transport of sensible 
heat in the vertical is a direct function of stability. 
Increased atmospheric instability ( 207) pol’) implies greater 
upward heat flow while conditions of strong stability 
(J6/) p</') will create downward heat flux. In the ordinary 


quasi-geostrophic two-level model with the static stability 
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Fig. 2.1. Vertical structure of the Lorenz (1960) two-level model 
of the atmosphere showing the isobaric information levels. A 
measure of the atmospheric static stability is given by 


(6. = SAGE = p,)- 


constant it is not possible to include the effect of a 
variable rate of vertical heat transport. As a result 
diabatic heating in such models is normally input at level 2 
(Figure 2.1) as a mean rate for the entire column of air 
extending from the surface to the top of the atmosphere. 
Since the whole atmosphere is essentially treated as one 
layer for heat input purposes the vertical diffusion terms 
can be ignored. Haltiner's (1967) vertical profile as given 
by (2.2.4) also precludes the need for vertical diffusion 


terms, although the profile is very dependent on the choice 
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of a value for the constant r which has no physical 


interpretation. 


2-3 Frictional Terms 

In Chapter 1 reference was made to two studies which 
reported investigations of the dynamic response changes 
created by inserting friction into a quasi-geostrophic two- 
level model. However, as used by Holopainen (1961), and 
Haltiner and Calverly (1965), the model embodied only 
surface-based frictional effects. In the current 
investigation those erictional effects created by small- 
scale turbulent eddies in the atmospheric flow well removed 
from the earth's surface will also be considered. 

The governing vorticity equation for the present model 
will be of a form similar to (2.1.14) except that a non-zero 
term will be inserted to account for frictional influences. 
Since it appears by way of the vorticity equation this 
frictional term will be treated as a rate of vorticity 
change caused by eddy momentum diffusion. If the term Fi 
represents the effect of friction in the model at an 
isobaric level j then, as with the approach employed for the 


sensible heating, it can be separated into two parts 
Fo = Neh (2.3.1) 
ein. “5 deere 


in which N and M are the turbulent diffusion terms 


representing the vertical and horizontal frictional effects, 


respectively. 


ohh Shadow t Beg 
OMe jcteen Oeoee Ay ive ‘le eee eee e 

g-bisihy (de nega 

02 be > dip ¥: ka 
7 ead cae ee! ue 

sa bg mm Poe si AR oe 

ays aes Aa ie «Ee goeeu Lat aaa H 

try jrnhre gp 278 ab 4 Sac te 

vel Lety, asa ole Angas 

Hig et i al Deri. 

nas 4ut kt au ‘we 

ss 

fauton (ent PERG AF EES S 408 erase og ‘pe Pua a sag 

EIey AV sISRBS hs ol eh ty. wae 4 oer am 

fb Peay Reet eae iomes ed Liat aed Ly ae 

(2 ape oer CI ert ty wATeS0EY Shoe ib ne 


ar: Hod é Linheeg ayes aa! ehoota’s 3p FIQRaS- ead 


od dokifaye isa. deeGerRA WAR) wtsW BE gana ay owen A } 


ac ake “Bee are roeperere a aan 7} series sid e 
a 
Ree ie , , : + ea) vA 4 As, ioe ’ i u Ss 4 A 
"Ua oS 3 s : a haa oan - 
° 2 re Eos \e ‘\, fi ah") Oy pts a - ; - = 6 


$714 cian lane At ihn ieee tan i i rad 


q 


By analogy with the thermal diffusion term R, the 


frictional effect due to vertical momentum diffusion will be 


given by 

Nea, d% 

ie oak: dp (25322) 
where ty represents the curl of the frictional stress acting 


on a horizontal surface and is defined as 
=> 
ma oV | 
x 
r- Vx ok (%) Bae, PELE 
2 Pe (2. 3. 3) 
for which K, is the vertical kinematic coefficient of eddy 
viscosity. Similarly, following closely the form given for 


the heating term S, the horizontal momentum diffusion tern 


will be written as 
IG x(GVY)-GV(VY) ae 


in which K, is the horizontal kinematic coefficient of eddy 


viscosity. 


2.4 The Model Equations 
Following the vertical structure as given in Figure 2.1 
the surface pressure value at the bottom of the current 
model is taken to be 100 kPa. The uppermost level is then 
given a pressure of 20 kPa as a representative value of the 
top of the troposphere. A common alternative is to let the 
upper level of the model correspond to the top of the entire 
atmosphere (p=0). However, for this study it will be 


assumed that the tropopause acts as a lid on the 
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tropospheric motions. The remaining information levels in 
the model are arranged such that each is separated by a 20 
kPa layer. Thus levels 1, 2 and 3 correspond to the 40, 60 
and 80 kPa frressure surfaces, respectively. 

The pare cnicd equations as applied at the different 


levels in the modelled troposphere are 


7(34)- #-0G8)-t0(F)-0 me 
3, T(Va)s V1- ie) ra (3) as (2.4.2) 
ai + Jb otf) 0% f- iShers Fees} 


where j=1 or 3. Equation (2.1.15) has been differentiated 
with respect to pressure and applied at level 2 to form 
(2.4.1). The diabatic thermodynamic equation and vorticity 
equation (friction included) as employed at levels 1 and 3 
are given by (2.4.2) and (2.4.3). The compatible form of 
the continuity equation remains as (2.1.12) with (201611) an 
inherent assumption. 

The equation of state for an ideal gas will also be 


used in the form 
pz prt (2.4.4) 


and (2.1.7) is understood to be an additional equation of 


state applicable to the model. 


The purpose of differentiating the linear balance 
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equation with respect to pressure is to derive a formulation 
which does not contain geopotential as a variable in the 
model. Replacement of oso using (2.1.11) and the 
subsequent use of (2.4.4) and (2.1.7) allows (2.4.1) to be 


rewritten as 
K-12 


2 7-Hf-(0%) fA) ans 
The tern i tote the change of the rotational 
component of the wind with changing pressure and can be 
thought of as a form of vertical shear. Therefore (2.4.5) 


is a form of thermal wind relation for the mid-troposphere. 
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CHAPTER 3 


NUMERICAL PROCEDURE 


If the earth's surface is assumed horizontal and the 
lowest level of the model remains fixed on the ground the 
proper lower boundary conditions are z=0 and w=dz/dt=0, 
where z is height. However, for convenience the conditions 


at the bottom cf the present model will be taken to be 
® =0 and P=p,= 100 kPa (3.1.1) 


These conditions state that the pressure at the bottom of 
the model remains fixed, but they do not imply a constant 
height field because z is not necessarily constant at the 
lower boundary. AS will be discussed in detail later in 
this chapter, level 4 of the model is not always found at 
the sea surface. However, the variation of sea-level 
pressure from 100 kPa will be.a perturbation quantity. 
Therefore, the lower boundary of the model will be treated 
as if it remained within a surface-based layer where the 
eddy stress term and the vertical heat flux are constant. 
No orographic variations are included in the lower boundary 
conditions. 


At the top of the model the condition that the 
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tropopause acts as a lid on tropospheric motions gives 
® =0 at p=p,=20 kPa (3e:1.2) 


To prevent the diffusion of heat and momentum across the 
tropopause the upper boundary conditions will also 


necessarily include 
h=0 and? =0 at p=20 kPa (ies) 


By adopting these boundary conditions it is possible to 
expand the friction and heating terms more explicitly. fThe 
expression for the vertical diffusion of sensible heat 
(2.1.6) is evaluated at level 1 using vertical finite- 


differencing to give 
A 4 
= pay» 
J Ap (3.1.4) 


where p represents a 40 kPa interval and h, is the vertical 


hee 
R aes 


heat flux at p=0 which by (3.1.3) is zero. Using (2.2.7) to 


replace h, (3.1.4) becomes 


24 LG 3-4, 
4S dp (4; HL SOLO 


The coefficient k,=p%g2K, /Ap? has been introduced into 
(3.1.5) for simplification. The same procedure can be 


applied at level 3 so that 


beh 
Sf 


eV 


(3.1.6) 


where the term h, is to be replaced as was done in (3.1.5). 


At the lower boundary of the model a constant surface heat 
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flux layer has been assumed. Surface heat flux will now be 


defined to be 
pee pglh- ie) A / | (3.1.7) 


for which C, is the surface drag coefficient. To convert 
(3.1.7) to a heating rate per unit mass in a form similar to 


that of (2.2.3) a constant C is introduced, where 


C 
Ge Ri ci PP E)) 


GP 


The resulting expression for R; is, therefore 
K, = “h Gp (3-4-4, )+ CVG-& ) (3.1.9) 


The entire approach as outlined above can be applied in 
a similar way to the vertical momentum diffusion terms. At 
level 1 the expressicn for the friction term in the vertical 


using (2.3.2) and (2.3.3) is 
2, -7, : %g , 
NV, =a, ae = k 7% -%) (3.1.10) 


where the coefficient k; =p2g2K, /Ape is analogous to k,, At 


the lower boundary the constant stress term is defined to be 


Apk, 2 
Gatien a / y (3212.19) 


for which k, is a constant. The term N; is found to be 
sh LVE-Y)-AI'Y 
WS AY . 3 AW 4 p3.1.12) 


It now remains to assign numerical values to the 


coefficients k k ; and K .° For the purpose of estimating 
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these values it will be convenient to assume K,=K,, and 
therefore k;=k,. In unstable cases there is some basis for 
the argument that the magnitude of K, is greater than K,, 
(Lumley and Panofsky,1964; Businger,1973). However, over a 
large range of Stability conditions K, and K,, have been 
found to be of similar magnitude. As well, the value chosen 
for K, (K,) in this study will be indicative of near-neutral 
stability conditions rather than of the unstable regime. 
Using J6yap2! ana values suggested by Deardorff (1967) K, is 
taken as being 6X10! m@ys. A representative 60 kPa winter 
temperature value at mid-latitudes would be between -18°C 
and -25°C based on mean January temperature data 
(Gates,1975). Applying (2.1.7) to determine density, 
possible values for k; and k, are then 0.9X10-® to 1.1X10-¢ 
s~i1, Estimates of k; as used in other two-level models 
range from 0.5X10-© s~! (Charney,1959; Kikuchi,1969) to 10-6 
srt) (Eggery1976)..% Forsthisiinvestigationyk;=k,=107% szivis 
used in all numerical comrutations. 

The most commonly used value for the coefficient k, is 
4X10-© s-! as employed by Phillips (1956) and Charney 
(1959). Derome and Wiin-Nielson (1971) have experinented 
with a surface friction coefficient which was altered 
between 4X107-® and 6X10-® s~!. The effect of increasing the 
surface friction coefficient on the behavior of their quasi- 
geostrophic steadyzstate model was relatively small. These 
investigators also studied the effect of having one constant 


surface friction coefficient over the land and a smaller 
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constant value over the oceans. The consequence on the 
model predictions was only of serious import when the ratio 
of the two coefficients had reached a value of 6 or more. 
At a ratio of 2 little difference was observed over assuming 
one uniform friction coefficient. The earth's surface will 
be considered to be uniform for the purposes of this study 
and the imposed surface fricticn coefficient will be taken 
to be 4X10-© s~! in most calculations. 

Other coefficients to be assigned numerical values are 
K, and K,. Once more, for simplification, it is presupposed 
that K=K In assuming Fickian diffusion for the 
dispersion of clouds and inert tracers in the lower 
atmosphere, the observational evidence as suggested by Hage 
(1964) and Bauer (1974) would indicate that the order of 
Magnitude of the horizontal diffusivity coefficient is 10S- 
106 m2/s for large-scale diffusion processes (i.e. those 
with a travel time of more than few days and an eddy space 
scale of greater than few hundred meters). Although the 
data refer to the horizontal dispersion of mass and may not 
be totally indicative of linear momentum or heat dispersion, 
they do provide a useful quantitative estimate of the order 
of magnitude which should be looked at. Phillips (1956) has 
employed a value of K,=105 m@/s in his modelling as does 
Kikuchi (1969). The value fixed for K, and K, in this study 
is also 105 m2/s. 

Surface sensible heating, as was indicated earlier, is 


dependent on the product of the constant C and the mean 
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surface wind speed V;. In order to assign a numerical value 
to C it is first necessary to make an estimate of the 
surface drag coefficient C,. Cressman (1960) has devised a 
general scheme for calculating surface drag coefficients. 
The coefficients are evaluated using two components, the 
first being a dimensionless constant of 1.296X10-3 for all 
Surfaces and the second is an additive quantity which is a 
function of the terrain roughness. Over water surfaces the 
total drag coefficient is resolved to be <1.3X107-3 while 
over land surfaces the variation in C, values is 1.3X1073 to 
8.0X10-3, with the larger end of this scale reflecting very 
mountainous terrain. If it is assumed that C,=1.3X10-3 for 
a surface wind of V,;=10 m/s (i.e. measured at anemometer 
height) then a numerical estimate of C as calculated fron 
(3.1.8) would be 8X10-3 to 9X10-3 J/(kg m °C). The product 
CV, then compares favourably with the value of 10-2 J/(kg s 
°c) used by Haltiner (1967) to represent the same quantity. 
As a matter of procedure the rate of sensible heat input 
will be varied when analyzing the model's dynamic stability 
characteristics. The heating rates used will remain within 


a range of 0 to 2x10-2 Jy(kg s °C). 


3.2 Perturbation Equations 


: 


A useful technique for representing conditions in the 
mid-troposphere is to create new variables which describe 
the vertical mean and differential quantities across levels 


1 and 3 of the model. For example, if b is any scalar 
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quantity then the averaged or mean value of b for the layer 
between the two levels is given by (b,+b,)/2. Similarly the 
differential value would be (b,-b,)/2. If this method is 
applied to (2.4.2) and (2.4.3) the equations take the 


following forms 


5 alle sales Clare ly a ig aamaanene ae 
2 He) aThya)«Ph-De (8) -(Y Haas 
Thy ef)i Th Th) ope 
Pedy O Vel) Sb) Vl Vt. Be B® aay 


Equations (3.2.1) and ae now represent the averaged 
thermodynamic and vorticity equations, respectively, while 
(3.2.2) and (3.2.4) are their difference counterparts. The 
variables as they now cccur are mage W726 Math Soe 
O= ( O+ PAV 2% c=(0- Q)/2 and 1 =(%,- Wize Vertical finite<« 
differencing has been used for those terms involving (2/dp) 
hae it has been assumed that at any given time 0 is 
representative of ene static stability over the entire 
troposphere (9 =Ojss +0). 

With the approximation that O = 9s 05) 72 (2.4.5) as 
applied at level 2 becomes 


ke 


ne V9 ” VE Vn) +0 (3.2.5) 
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while the continuity equation for the model using the upper 


and lower boundary conditions on &,can be expressed as 


| )-(32)|. £ ny7eel (3.2.6) 


In order to analyze the dynamics of this system a 
perturbation approach is used to linearize the equations. 
The perturbation method assumes that atmospheric motion 
consists of small fluctuations superimposed on a basic or 
steady-state flow. The fundamental premises are (Haltiner 
and Martin, 1957): 


(a) the steady-state motion satisfies the 
governing equations; 

(b) the overall] motion (steady-state plus 
perturbation) also satisfies the equations; 
and 

(Cc) Squares (or higher orders) and products of 
perturbation variables appearing in the 
equations may be neglected in comparision with 
first crder terms. 


For instance 


Y. V, yp’ | | (3.2.7) 


preven denotes the total value of the quantity, V is the 
steady-state value and } refers to the perturbation value. 
The same representation can be applied to (ee Oo and Me 
The steady-state flow as is suggested by its name is 
not a function of time. Further, it will be assumed that 
such flow is zonal and non-divergent in behavior. Thus, 


both % ana & are necessarily zero. The basic state coraae 
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6,7 ana O are to be regarded as linear functions of y alone 
while the perturbations are to be treated as simple 
Functions, of x, y and. t. 

For sensible heat exchange the earth's surface is 
regarded as on infinite heat source (or sink) implying that 
its temperature field remains constant and that T, tends to 
adopt to Tye Therefore, under steady-state conditions the 
surface air can be assumed to be in thermal equilibrium with 
the underlying surface. As a result T=, and the heating 


rate is dependent on surface air temperature fluctuations 


Oban) k, (-7’) (3.2.8) 


Variations in V, may be neglected since T, is already a 
perturbation quantity and so k, =CV, « 

If the potential temperature at the lower boundary of 
the model is linearly extrapolated from the 60 kPa level as 
O=6-20, then by assuming a sea-level pressure of. 100 kPa 
the atmosphere-ground temperature difference is 1, = 26-20. 
This assumption is not totally valid, however, and a 
correction term is Peanieed to describe the possible 
variations in the Beaticvel pressure from 100 kPa. 

An estimate of the pressure deviations at z=0 can be 
obtained by using the general form of (2.4.1) and 
substituting for the geopotential with ?=gz. For a ticst 
approximation the term Ve VOY dp) is neglected and the 


remaining partial derivatives are replaced by finite 
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numerical estimates giving 


sz £ SF 

bp Ber Sp ANE PAY 
where oz and 9) represent changes in the height and 
streamfunction values over some very small pressure interval 
Sp. Perturbations in the lower boundary boundary stream 
function values are given by Y, V2. Hence, assuming syi¥ 
and using (3.2.10) the variation in the 100 kPa level height 
is Sz=t(P-21)/. Moreover, from the hydrostatic assumption 
(2.1.11), the change in pressure in going from 100 kPa to 
sea-level pressure is 5 p=-pt (V=27) Figure 3.1 illustrates 
a Situation where Y-21>0 and the sea-level pressure is 
greater than 100 kPa. 

From (2.1.7), (201.211) and (2.4.4) it is possible to 


derive the following expression for the atmospheric lapse 


oT [20 a\/p \x 
dz ; (38 4 (£) "| Wp (3.2.10) 


in which Y=9/c, is the dry adiabatic lapse rate. Replacing 


rate 


the partial derivative dT/dz with its finite representation 
$a/0z and using d0/bp=-20/4 p for the modelled stability 


parameter gives “st 
Sy eee ly, r ne (EY (3.2.11) 


where 5T is now the temperature decrease in going from z=0 
to the height of the model's lower boundary 6z. The total 


variation in the surface temperature is 
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Fig. 3.1. Deviation of the sea-level pressure field from 100 kPa. 
Sea-level pressure values are estimated using 100 - ép (kPa), 
where Sp = - pgéz = - pf(p - 2n). The vertical cross-section 
shown displays a situation where the 100 kPa level is found above 


sea surface. 
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To a very good approximation ((100- Sp) 7100) ~ 1 and as a 
result may be dropped from (3.2.13). The total surface 


heating rate then takes the perturbation form 


-k [a tlhe 6° gh an Nee Zp 2) (3220113) 


In addition to having no surface heating under steady- 


state conditions, it will also be presupposed that there is 
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no sensible heating from diffusion in the basic flow. This 
has already been accomplished for the horizontal diffusion 
term by assuming that the steady=state quantities are linear 
on the y-axis and constant along the x-axis. In the 
vertical diffusion term it is achieved by setting the 
countergradient s6/dp equal to -20/Ap in the model. 

Vertical diffusion at level 1 as expressed in (3.1.5) is 


then 


26 
is - kc, Ae Ap A | is “ake or. (3.2.14) 


6 P 


A similar substitution can be made in (3.1.9). The value of 
O used for the modelled atmosphere is based on January 
average temperature fields at 40 kPa and 80 kPa for the 
Northern Hemisphere (Gates, (1975)). Over the hemisphere 
and, in particular, over the North Pacific the mean value of 
the static stability is approximately 22 to 23°C for the 40 
kPa layer. This would give a countergradient aa)z of 4 to 
5°C/km. Typically, estimates of [in the atmosphere vary 
from 0.65°Cykm under clear skies to roughly 5°C/km under 
cloudy conditicns (Deardorff,1967). 

The linearized versions of (3.2.1)-(3.2.5), after 


substituting for w, using (3.2.6), take the form 
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(3.2. 16) 


(3.2.17) 


(3.2.18) 


(3.2.19) 


where U=-dY/dy represents the mean zonal wind and U=-7/0¥ 


is the vertical wind shear of the basic state. 


coefficients a, 


= (4/5) ~k- (3/5) ~*=-0.233 and a, =(4/5) —“/2=0.533. 


3.3 Stability 
In order 


perturbations 


Analysis 


to investigate the stability of those 


The 


(i=1,2,3) are; a,=(2/5)—*+ (4/5) —k=2.635, 


in the baroclinic atmosphere depicted by the 
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model, harmonic-wave soluticns of the following form are 


assuned 


ve Vyas lene ef | H> 7 (y)exo iB (x ac ) 


O’- Gl.) oxp he ) ae © (x)exp iP (,-ct)| (3.3.1) 


X- X lye] inl) 


where m= 27/1 y is the wavenumber and L, is the wavelength on 
the x-axis, c is the phase velocity of the wave, Vi= -1, and 
Yon OO. © PGi are the wave amplitudes. 
The latitudinal dependence of the perturbations enters 
through the amplitudes which will be given a simple 


meridional variation 


Vy : Vos oy) eee 


where nS2 1/1, is the wavenumber and Ly is the wavelength in 
the y-direction. Substitution of these solutions into 
(3.2.16) -(3.2.20) leads to a set of homogeneous equations 


which can be represented in matrix form as 
(BJ{ A ]=0 (3.3.3) 


Here {A] is the amplitude matrix 


WY Ds, 
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and {B] is the coefficient matrix defined as 
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In the coefficients G= ( Yj-2ego/ p) and 4=df/dy is the 
meridional derivative of the Coriolis parameter. Use has 


also been made of the following relations 


26. -s*fu 
dy ce uk (3. 3.4) 


a K 
5 (1-x) ae (3.3.5) 
which were found by assuming a constant U,; in the mid- 
troposphere and then requiring the zonal flow to satisfy 
(3.2.20) « 

A necessary and sufficient condition that (3.3.3) has a 
non-trivial solution is that the determinant of ({B] be zero. 
Expansion of the determinant and grouping of the terms in 
decreasing orders of (c-U) gives a cubic frequency equation 
which can be solved for a variety of wavenumbers, surface 
heating rates and steady-state quantities. 

The phase velocities which are obtained in the 
solutions to this cubic frequency equation are typically 
complex values, c=c,+ic., in which c, represents the real 
portion and c_ the imaginary portion of the wave phase 


velocity. A perturbation wave is, therefore, of the form 


Vs ¥ cos(ny) exp (met ) exp (iim (x -c,t)) (3.3.6) 


Such a waveform propagates in the west-east direction at a 
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speed C, with a time-dependent numerical amplification (or 
damping) factor of exp(mc_t). More specifically it is the 
value of c, which determines the stability condition of the 


wave. The stability classifications are given in Table 1. 


Table 1. Stability condi tions for the harmonic wave- 


forms as a function of the imaginary portion of the phase 


velocity. 
STABILITY CLASSIFICATION 
UNSTABLE : Cy <O , amplifying wave with time 
NEUTRAL: Cy =0 , waveform unchanging 
STABLE : cy >0-., damped wave with time 


The rate of amplification of an unstable wave is 
dependent on the exponent mc,. A useful way of representing 
this growth rate is by a quantity known as the e-folding 


time (t,) which is defined by 


ae Z 
tm, * 296 (3.3.7) 


The value of t. is the time required for an unstable wave to 
increase in amplitude by a factor of e=2.71828. The 
implication of a small e-folding time is that the waveform 


is very unstable and amplifying rapidly. 
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3.4 Wave 


Dynamics 


=P qe ee 


When studying the phase and amplitude characteristics 


of the perturbations it 
solutions in a somewhat 
each value of the phase 


the frequency equation, 


by (3.3.3) can be solved for values of ie Age on 


is useful to express the harmonic 


different form. Corresponding to 
velocity c, (k=1,2,3) as given by 
the system of equations represented 


o and 7, 


if initial conditions for each perturbation wave are 


specified. 


The procedure is then to write the general 


perturbation waveforms as sums of components 


kz 


bey 


ot 


cos (ny) cxp im(x-¢t) 


ne >, i, ces (ny) exp im (x-6t) (3.4.1) 
O- 2. G, cos ny) exp im (x-¢£) 

Oo’ va cos (ny) ox trea (x af) 

UT migpepmtt) 


It is convenient to substitute for the terms 1, and x, 


from (3.3.3) using 
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where [E]= vy, , and ({DJ= dt a? 0 
4, be dg a 4a§ 
0, de edi de 
The coefficients a) (j=1,2,+++,8 k=1,2,3) are 
a= (¢-U/) ob a 
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The initial conditicns specified at t=0 in the 
solutions to (3.4.4) are of the type 
/ A 
V [£-0) = iv COS (74) exp (ix) 
/ A ; 
a) (t-0) = 0, COS (rq) exp (im Xar a) (3.4.5) 
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in which ; @. and 05 are the initial amplitude values and 
ae and € are the initial phase angles of Bana © with 
respect age For example, £, <0 implies that the potential 
temperature wave initially lags the stream function, while 
é, >0 means that the potential temperature wave leads. If 
E,=0 the two waves are initially in phase. 


It follows from (3.4.1) and (3.4.5) that 
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ay, G, ; Q exp (té,) (3.4. 6) 


Furthermore, from (3.4.4) 0, and OF can be expressed in 
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terms of } as 


het kot 
3 3 A (3.4.7) 
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Knowing the initial amplitude and phase values, the 


phase velocity Cc, and parameters such as U, f, er and 


A 


others it is possible to determine Me and then, in 


A A “ Pa aN 
succession, the values of ae fs ; GO, y 6, and w, with the 


component amplitudes for the vertical velocity determined 
from using the model's continuity equation and the velocity 
potential values x. It is also useful to derive wave 
components representing surface sensible heating ch, ) from 
the existing waveforms by employing (3.2.14). 

Although each perturbation is written as a sum of three 
components, a single wave representation can be made which 
has its amplitude and phase angle as functions of time. 

From (3.4.1) and by writing the phase velocities and 
amplitudes in their complex form the mean perturbation 


stream function may be expressed as 


ion (Pit) cxtrydep| met + im (rg) (3.4.9) 


In complex polar coordinates this is equivalent to 


met <éla-met) nx 
"6 (t ‘ee (3.4.10) 
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in which i? | = feny 2)1/2 Eeprezents the modulus of the 
complex vector pap and q -arctant  /f 1 is the vector 
angle from the positive real axis. SETI EE) the mnodulus 
and amplification factor into one term H+ /P] expinc,t 
and defining q ty were allows (3.4.10) to be rewritten as 


a single wave solution of the type 
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where the amplitude uy ist 


A kat R kot 
<— y (3.4.12) 
Soe Ay ay * ali z 
29 Cos 4 ae O ZY, sing. 
ks2 
and the phase angle \) is determined from 
2584 Ps 
W Ca) 
= arctan ker k Zz 
( 2 y* é ) 
k & ‘ | 
% § (3.4.13) 
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In the same manner the remaining perturbation 
/ / 7, if, 

quantities Boys Oligg Gee u and “~ may be described in a form 
egupie to that of (3.4.12) with their respective amplitudes 
A Ve 4 A Q 

Q, , 17, and phase angles G.o 4 ¢ Np ¢ 7p 
Calculations were carried out to determine these values as 
functions of time and for varying initial conditions as well 


as for different basic flow properties. 
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CHAPTER 4 


NUMERICAL RESULTS AND DISCUSSION 


4.1 Preliminary Comments 

The present model, like that of Haltiner (1967) 
incorporates a surface sensible heat flux but, unlike the 
latter, it also allows for the inclusion of surface friction 
effects, heat and momentum diffusion terms and a simple 
meridional variation in the perturbation field. The 
numerical investigations of the dynamic response 
Characteristics of the model for different states of current 
flow are discussed first. These flow regimes are 
indentified in the following manner: 


(i) Model 1. Adiabatic flow with no frictional 
vorticity terms included. 


(ii) Medel 2. Sensible heating (including the 
redistribution of heat by thermal diffusivity) 
with no frictional vorticity terms. 


(iii) Model 3. Adiabatic flow with frictional 
vorticity terms included. 


(iv) Model 4. Sensible heating with frictional 
vorticity terms included. 


As an initial simplification it is assumed in each of 
the four ‘models' that the perturbation field is nota 
function of latitude. Later, a latitudinal dependence is 


incorporated in conjunction with model 4 and the 
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calculations are redone. This particular version of the 
modelling is simply referred to as the ‘full model’. 

The latter part of this chapter focusses on the 
behavior of the longest waves. Of specific interest are 
those features of long waves which can be applied to 


blocking ridge development in the North Pacific. 


In order to compare the various model runs the steady- 
state atmospheric quantities were held constant unless 
otherwise noted with U,=10 m/s, f (45N) =1.0312xX10-* s-!, 
AU45N) =1.61X10-12 m-ts~! and O=12°C. The phase velocity 
solutions of models 1 through 4 were then obtained from the 
corresponding cubic frequency equation. Since the frequency 
equation is solved in terms of (c ,~U) it should be kept in 
mind that the real porticn of any solution represents the 
phase velocity of the perturbation wave relative to the 
speed of the mean westerly wind. 

Figures 4.1a and 4.1b show both the imaginary and real 
portions of the solutions (c,-0) as a function of wavelength 
L, for nodels 1 and 2, respectively. Each case exhibits a 
neutral or neaE-neutral solution (c,~0) arbitrarily denoted 
Ii, a damped solution (c.<0) designated curve I2 and an 
unstable solution I3 (c, 20). The corresponding real parts 
to these curves are shown as R1, R2 and R3. 


Features apparent in the adiabatic-frictionless 
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Table 2. Real (R3) and imaginary (13) portions of the phase 
velocities over the unstable waveband of Fig. 4.la and the 


corresponding e-folding times (t.). 


Wavelength (km) R3 (m/s) 13 (m/s) te (days) 

1000 =7 220 0.00 

1500 moe 2.05 14:35 
2000 -2.80 4.50 0.82 
2500 “25/5 5699 0.77 
3000 =3.00 6.95 0.80 
3500 -3.49 7253 0.86 
4000 -4.14 7.82 0.94 
4500 -4 93 7.86 1.06 
5000 -5 83 7.66 1.20 
5500 -6.86 7.16 1.42 
6000 hoo 6.27 1.76 
6500 -9.23 47h Pi ctov2 
7000 =1i.05 0.00 


solutions of Fig. 4.1a are that the unstable and damped 
waves propagated at the same speeds over the entire range of 
unstable wavelengths. Furthermore, these speeds were always 
slower than that of the mean westerly wind speed. The 
neutral wave in turn always travelled faster than the zonal 
wind speed as indicated by the R1 curve being positive at 
all wavelengths. 

Introduction of sensible heating to the model as shown 


in Fig. 4.ib has diminished the amplitude of the I3 curve at 
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Table 3. Real (R3) and imaginary (13) portions of the phase 
velocities over the unstable waveband of Fig. 4.1b and the 


corresponding e-folding times (t.) 


Wavelength (km) R3 (m/s) I3 (m/s) t. (days) 
1000 = Je 27 0.06 30.70 
1500 -4.59 Jails 2.45 
2000 -3.07 2.81 ey 
2500 =2855 4.20 1.10 
3000 -2.72 5.07 1.09 
3500 ~3.27 5.47 1.18 
4000 -4.09 5.50 1.34 
4500 -5.18 Gea? 1.59 
5000 -6.61 4.67 1.97 
5500 -8.47 3.93 Pe 
6000 -10.82 ely, 3.49 
6500 12657, 2.54 4.7] 
7000 -16.60 2.06 6.26 
7500 | -19.85 1.70 8.13 
8000 273°28 1.42 10.38 
8500 -26.90 1221 12.94 
9000 -30.69 1.05 15.79 
9500 -34 66 0.91 19.23 

10000 -38.82 0.81 22.74 
10500 43.15 0.72 26 . 86 
11000 -47.68 0.64 31.66 
11500 -52.39 0.58 36.52 
12000 -57..30 0.53 4y.71 
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lower and intermediate wavelengths, but at the same time has 
extended the band of unstable wavelengths to beyond 10,000 
km. Other changes from the values cf model 1 are that the 
I2 curve now shows much stronger damping at longer 
wavelengths and the R2 and R3 curves no longer coincide over 
any interval. The coefficient for surface heating in this 
case was chosen as k,=1072 J/(kg s °C) and, unless otherwise 
noted, this value was used in all calculations where heating 
was incorporated. 

Haltiner's (1967) findings using the same value of 
surface heating and a vertical wind shear of 10 m/s over a 
25 kPa layer show qualitatively the same shape as the curves 
of Fig. 4.1b. Tables 2 and 3 give the vaiues for the 
unstable waveband at every 500 km for models 1 and 2 from 
the present study, while Haltiner‘'s results are found in 
Tabie 4. The latter two tables are quite comparable over 
wavelengths from 3500 to 8000 km. Outside of this range the 
present model gave a slightly more unstable wave as 
indicated by the c, values and the computed e-folding times. 

The results of incorporating fricticn are presented in 
Fig. 4.2a. It is apparent that there are two regions at the 
extremes of the unstable waveband where increased 
instability occurred when compared with that of the 
adiabatic-frictionless case. The amount of broadening of 
the waveband at longer wavelengths was very small, however, 
when compared with Fig. 4.1b. It is aiso evident by 


comparison of Tables 2 and 5 that, except within these two 
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Table 4. Real (R3) and imaginary (13) portions of the phase 
velocities over the unstable waveband of Haltiner (1967) for 
sensible heating with no friction. (k, = 10 2 J/(kg 5 cys 


r = 2, U. = 10 m/s ando = 15°C over a 25 kPa layer) 


T 
Wavelength (km) R3 (m/s) 13 (m/s) i (days) 
1000 -8.39 0.00 
1500 -6.3] 0.24 Lies 
2000 -4 45 2.06 lew 
2500 -4 03 3.49 lez 
3000 -4 06 4A 1.24 
3500 -4 uh 5.03 1828 
4000 <5. 12 5-30 139 
4500 -6.07 5.50 1.56 
5000 oy ops 18) 5.04 1.83 
5500 -8.85 455 2.28 
6000 -10.80 3.88 2.85 
6500 = 3719 3:10 3.86 
7000 -16.04 Z5ih 5 AA 
7500 =A O27 1.78 7.76 
8000 -22.76 1.36 10.84 
8500 -26 45 1.06 14.77 
9000 =30).511 0.85 19.50 
9500 -34.34 0.69 25.36 
10000 -38.54 0.58 31.76 
10500 ~42.92 0.49 39.47 
11000 -47 47 0.42 48.24 
11500 = 52521 0.36 58.84 
12000 -57.14 0.31 Ate 
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harrow regions, the amplitude of the I3 curve has everywhere 


been reduced. 


Table 5. Real (R3) and imaginary (13) portions of the phase 


velocities over the unstable waveband of Fig. 4.2a and the 


corresponding e-folding times (t,). 


Wavelength (km) 


1000 
1500 
2000 
2500 
3000 
3500 
4000 
4500 
5000 
5500 
6000 
6500 
7000 
7500 
8000 


R3 (m/s) 


60 
66 
102 
-03 
-50 
ay 
ahs 
SoZ 
95 
18 
54 
Bus) 
.66 
a0 


I3 (m/s) t. (days) 
-8.23 
1.45 1.9] 
3.68 1.00 
5.11 0.90 
6.04 0-92 
6.59 0.98 
6.86 1-07 
6.88 2 
6.65 1.39 
6.15 1.65 
5.29 2.09 
3293 3.05 
Lo 6.75 
0.28 49.34 
<0 


—————————— —_______ _______________________ EEE 


The response curves of model 4 containing both sensible 


heating and friction are given in Fig. 4.2b. A longwave 


cutoff is now found near 8500 km, roughly 1500 km longer 


than the cutoff shown for adiabatic-frictionless conditions. 


It is also significant that over the unstable waveband the 
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R3 values of Fig. 4.2b imply much lower propagation rates 
than found initially in Fig. 4.1a. For instance, a mean 
wind of U=15 m/s would mean that every unstable wave in 
model 1 would propagate eastward as it amplified. With the 


inclusion of both sensible heating and friction this value 


Table 6. Real (R3) and imaginary (13) portions of the phase 
velocities over the unstable wave band of Fig. 4.2b and the 


corresponding e-folding times (t,). 


Wavelength (km) R3 (m/s) 13 (m/s) t. (days) 

1500 =36 bil <0 
2000 -2.81 UATAw De Ah: 
2500 =e 3.08 1.50 
3000 -3.89 3.94 1.40 
3500 -4.89 4.3] . 1.50 
4000 -6.14 4, 33 1.70 
4500 =/e7 1 4.06 2.04 

- 5000 | -9.63 3.56 2.59 
5500 cal fl Pets 2293 3.46 
6000 -14.53 2.25 4.91 
6500 “17 oe 1.61 7.44 
7000 -20.56 1.04 12.40 
7500 =237.91 0.54 ja ESS) 
8000 -27.45 0.10 147.37 
8500 31319 <0 


for the mean wind would result in a stationary unstable wave 


of about 6000 km (Table 6) while at longer wavelengths 
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unstable wave growth could cccur as the wave retrogressed 
westwards. The implications of this will be investigated 
further in the discussion on blocking ridge development. 

The rate of surface sensible heat exchange and the 
degree of surface friction have as yet been held constant. 
Figure 4.3 illustrates the change in the unstable curve for 
model 2 under various rates of surface heating. The basic 
flow parameters were unchanged from the values used earlier. 
It is seen that increasing the surface heat input resulted 
in a further reduction of the curve amplitude at lower and 
intermediate wavelengths. Conversely the region at longer 
wavelengths exhibited increased instability. It is also 
noticeable that the region of enhanced instability between 
successive increments in the value of k, was gradually 
shifted to longer wavelengths, indicating that further 
increases in the heating rate would resuit in greater 
instability at extremely long waves only. 

Increasing the degree of surface friction through the 
coefficient k, in model 3 resulted in only marginal 
increases in the amount of instability at the edges of the 
unstable waveband. However, in changing k,; from 4X10-© to 
10-S s-1 a fairly significant reduction in the amplitude of 


the 13 curve did occur at intermediate values of the 


wavelength range. 


4.3 Variation of the Steady-state Flow Parameters 


The response values of the perturbation field are 
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affected if certain properties of the steady-state flow 
Change. This is easily seen in the case of a change in the 
zonal wind U. Since U appears as part of the real portion 
of the frequency equation solution it does not influence the 
Stability of the waves. Its main influence is to define the 
*absolute' speed of the waves on the x-axis and a change in 
U results in a corresponding shift of these values. This is 
important because the value assigned to U will play a major 
role in determining whether it is possible to have a 
stationary wave within the unstable band of wavelengths. 

The mean vertical wind shear U, is a measure of the 
baroclinicity of the basic state and the numerical results 
indicated that the greater the shear the greater the 
possible instablity of the perturbations in the flow. The 
solid curves of Fig. 4.4 represent the change in stability 
experienced in model 4 under the same conditions as employed 
earlier except that U; was varied, With greater wind shear 
the instability of the I3 curve was increased at all 
wavelengths. The effect on the real portion of the unstable 
solution was to decrease the speed of the waves by a few 
meters per second. In general very little instability was 
observed in any of the model types at wavelengths greater 
than 6000 km until the wind shear had exceeded 5-6 m/s. 
Thus, at low mid-tropospheric wind shear vaiues the 
possibility of longwave instability is small. 

Figure 4.4 also shows the effect on the I3 curve of 


changing the mean static stability. As the mean stability 
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of the atmosphere was increased the instability of the 
perturbation field was reduced. This reduction is much more 
Significant at the lower end of the wavelength range and 
became marginal at larger L, values. The relative 
propagation ate of shorter waves was also found to decrease 
as 0 was increased, but once again this effect diminished at 
longer wavelengths to the point where the wave speeds were 
virtually unchanged. 

Although the perturbations have been assumed to be 
latitudinally independent, variations in latitude may still 
influence the numerical results through the Coriolis 
parameter and its meridional derivative J. Therefore, the 
values of these two quantities were varied so as to 
represent different latitudinal x-axes. At lower latitudes, 
as suggested by a comparison of the unstable curves in 
Fig. 4.5, the perturbation field was found to be 
considerably mcre stable and the range of unstable 
wavelengths was reduced significantly. In addition, 

Fig. 4.5 shows that as latitude decreased the perturbation 
waves had a greater rate of retrogressicn with respect to 
the mean zonal Siete 

4.4 Perturbation Field as a Function of Latitude 

At this time the previous assumption that the 
perturbation field is independent of latitude is relaxed and 
the simple cosine variation of (3.3.2) is employed. Unlike 


wavelengths on the x-axis which have been expressed in 
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Fig. 4.5. Effect of latitude variations on the stability and 


propagation rates of unstable waves. The stability curves C73 and 


relative wavespeed curves CR3 - U were computed using model 4 with 


UL = 15 m/s. Latitude is distinguished by the bracketed values 


assigned to each curve. 


s~1, g(60N) = 1.14 x 1021 m!s 1, £(30N) = 7.292 x 105s } 


In the calculations f(60N) = 1.263 x 10-4 


and 6(30N) = 1.98 x 10 44 mis}, 
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distances of meters or kilometers a meridional wavelength 
will be given in terms of degrees of latitude. In all cases 
the origin in the north-south direction will be taken as 
45N. The interpretation of Ly=120° is) that: of a 

peetar bation eranediine in an east-west direction having its 
maximum absolute amplitude at y(45N)=0 and an amplitude 
which approaches zero at latitudes 75N and 15N. 

The computations were redone using the same numerical 
values as utilized earlier for the steady-state quantities. 
The full model was run for varying values of n=27/Ly « 

Figure 4.6 shows the phase-velocity response curves (c, -U) 
at latitude 45N when Ly =90°. The general shapes of the 
curves were Similar to those shown in Figures 4.1 and 4.2. 
A comparison of values in Tables 6 and 7 indicates that 
unstable long waves have greater rates of growth and travel 
at slower propagation rates after the 90° latitudinal 
dependence has been introduced into the modelling. 

Figure 4.7 illustrates the resulting unstable solution 
for a range of Ly values between 60° and oo (L, =e implies no 
latitude dependence). The changes that occurred in the c_, 
values in response to varying values of Ly were similar and, 
therefore, not all of these curves are shown. In going from 
Lh aN se OS the amplitude of the c,, curve showed a 
Marginal decrease in amplitude up to a wavelength of L,=5000 
km. At wavelengths L 25000 km instability increased and, as 
a result, the unstable waveband was extended to include 


longer waves. The relative travelling speeds of the 
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Fig. 4.6. Phase velocity response for the full model incorporating 
a 90° latitudinal dependence (Ly = Ga)" Curves RoR, and R, are 
the real parts and I, .1, and I, the corresponding imaginary parts 


to the frequency equation solutions (cy) =o)» kom il s2.3: 
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Table 7. Real (R3) and imaginary (13) portions of the phase 
velocities over the unstable waveband of Fig. 4.6 and the 


corresponding e-folding times (t.). 


Wavelength (km) R3 (m/s) I3 (m/s) a (days) 
1500 -4 34 <0 
2000 -3.18 1.54 239 
2500 -2.80 2.84 1.62 
3000 =3.03 3.68 1.50 
3500 -3.58 4.10 ie 7 
4000 -4 3) 4.20 | PAS A 
4500 Uys 4.08 2.03 
5000 -6.14 3.8] 2eaZ 
5500 Cay fen a 3.43 2.95 
6000 -8.39 299 5.70 
6500 -9.64 2.54 4.7) 
7000 -10.94 Pell Geist 
7500 = 12.25 P70 8.13 
8000 =13255 1233 11.08 
8500 -14.82 1.00 15.66 
9000 -16.05 0.69 24.03 
10000 -17.24 0.41 hh 93 
10500 -18.37 Cris 128.94 
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Fig. 4.7. Effect of variations in the meridional wavelength L on 
the stability and propagation of unstable waves. Numerical results for 


the stability curves c_, were computed using the full model and 


I3 
L = 60°, 90° and ~ as indicated by the bracketed values assigned to 
each curve. For the relative wavespeed curves the scope of YY values 
was expanded to include Hy = 120° and 150n- The curves for “8, = 0 


are interpreted as having no latitudinal dependence and are, there- 


fore, analogous to curves of model 4. 
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unstable waves were affected in a similar manner. Up to 
values of L, =4000 km the propagation rates of the waves were 
almost unchanged. At values greater than 4000 km reductions 


in the magnitude of L, resulted in the waves experiencing 


y, 
propagation rates which more closely approached the speed of 
the steady-state current. 

If y is taken to be non-zero in the model then the 
absolute amplitude of the perturbation field is reduced. 
This in Gesele has little bearing on the dynamic stability 
of the flow if the ratio y/ly is less than 1/4 (i.e. if 
y/L,~ 1/4 then tan(ny)>~© in the coefficient b!%). A value 
of y#0, however, implies higher or lower latitudes and the 
Main changes in the relative phase velocity curves occur 


because of changes in f and 8 aS waS seen previously in 


Fig. 4.5. 


4.5 Amplitude and Phase Relationships 

Phase velocity solutions to the frequency equation 
determine the stability features of the perturbation field. 
However, stability curves do not show directly the 
interactions of the three wave solutions or the 
relationships between the various perturbation quantities. 
The initial-value technique outlined in Chapter 3 will now 
be used to extend the analysis to include such interactions 
and interelationships. The following parameters were held 
constant: U=20 m/s, U, =10 n/S, L,=90°, kAsi0s? J/(kg s °C), 


=129C, k,=4X10-6 s-t, K/=K,=105 m2/s and k.=k =10-6 s-1 
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With ft and 4 evaluated at a latitude of 45N. In addition, 
the initial amplitudes assigned to the stream function, 
potential temperature and static stability perturbation 
waves were held constant at =9.806X10° n2/s, 6 =10°C and 
6, =2°C. In order to obtain estimates of the magnitude of 
the stream function it is useful to assume from the linear 
Re pA 
balance equation that ~=biyg where the mean Coriolis 
parameter can be taken as being approximately 107* s~1. 
Therefore, the averaged value of the stream function can be 
though* of in terms of a geopotential height with an initial 
amplitude of 100 m. In the future y will be expressed in 
its geopotential ‘equivalent in units of meters. 

Figures 4.8 and 4.9 show the characteristic curves for 
a wavelength of 4000 km taken from model 4. The initial 
phase angles €, and € in Fig. 4.8 were both taken to be 909 
meaning that the potential temperature and stability waves 
lagged behind the stream function. These angles were 
changed to -90° in Fig. 4.9. The 4000 km wavelength is 
indicative of the middle range of the instability waveband 
where growth rates are quite large. 

The rate of amplification for the different 
perturbation quantities as displayed in Fig. 4.8a was quite 
rapid during a 72-hour period. It is also apparent that, 
except for a time during the first 12 hours when the 
amplitude of the static stability briefly decreased, the 


amplification rate of each perturbation increased in time. 


Figure 4.8b displays all phase angles in degrees. If a 
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Fig. 4.8. Perturbation amplitudes and phase angles as functions of 

time for model 4 with L_ = 4000 km, U = 20 m/s, U; = 10 m/s, cy = 90° 
ree) : : ‘i l 70 Saye) 

and €. = 90°. Amplitudes in (a) are v, (10 m) , 6, ( C)s oA ( C)s 

wa(10 * kPa/s) and h,(10 ? J/(kg s)). Phase angles in (b) are in 


degrees. 
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curve slopes towards lower (negative) phase angle values the 
wave in question is progressing in an easterly direction. A 
Slope to higher (positive) phase angles thus popeesents 
retrogression westwards. All of the perturbation waves 
depicted in nie figure travelled eastwards and the short 
interval when om is damped is seen to correspond to a time 
when the static stability wave advanced toward the strean 
function. 

The growth curves for Fig. 4.9a are slightly changed 
from those cf the previous case. The stream function was 
damped initially for 24 hours while the potential 
temperature wave advanced to within a lag of approximately 
130° (Fig. 4.9b), after which the stream fuction began to 
amplify. In addition, the former wave's eastward motion was 
reduced substantially between 12 and 24 hours until the 
potential temperature wave had moved into this appropriate 
lag position. The potential temperature wave can be 
interpreted as a ‘thermal' wave since it uniquely determines 
the value of ny which is representative of a geopotential 
thickness in the same manner that YW" is representative of 
height. Thus, the chewed wave necessarily lags the 
geopotental height wave in order for amplification of the 
height field to occur. Static stability increased during 
this initial time period and then subsequently showed a 
short period of damping. A large drop in the surface 
heating amplitude was potiesd during the first few hours as 


this wave moved rapidly from a trailing position behind the 
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Fig. 4.9. Perturbation amplitudes and phase angles as functions of 
time for model 4 with L. = 4000 km, U = 20 m/s, Uy = 10 m/s, ee 90 
ra KO : rd A 1 n (0 “2 #6 

and €, = 90°. Amplitudes in, (a) are (10 m) , 6, ( a oA ( ae 


wy (10 kPa/s) and h,(10 ? J/(kg s)). Phase angles in (b) are in 


degrees. 
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stream fuction wave into a phase lead. 

Many runs were made with L,=4000 km for different 
values of Ete and with various model types. Bie. 
adiabatic-frictionless results of model 1 corresponding to 
the two cases already discussed have qualitatively the same 
growth and phase relations as when friction and sensible 
heating are included. In general, however, model 1 gave the 
greatest amplification of the waves during the 72-hour 
period. 

Certain phase difference values showed characteristic 
changes depending on whether heating or frictional effects 
were included in the modelling. As an example, in Fig. 4.8b 
the difference eal, at the end of 72 hours was 56°. In 
model 2, using the same initial conditions, this phase 
difference was 63° while model 1 gave a final value of 739. 
Thus, friction and sensible heating had the effect of 
allowing the thermal wave to approach a position somewhat 
closer behind the geopotential height wave during 
amplification. Other runs using different initial phase 
conditions values also supported these trends between the 
various model types. 

More pronounced changes were seen in the values of 
3 


stream function in 12 hours after which it slowly moved to a 


vt, - In Fig. 4.8b the static stability overtook the 


56° lead by t=72 hr.. Adding a 90° latitudinal dependence in 
the full model did not alter this value and dropping 


friction from the modelling resulted in only a 5° decrease 
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in the phase difference. Model 1, however, had a final 
phase lead value of 114°, a rather significant change which 
can be attributed to the lack of sensible Relinen | Once 
more this trend was supported in other runs with various 
initial Ponas tiene: 

The surface heating and vertical velocity perturbation 
fields showed evidence of having consistent phase relations 
with the stream function. Surface sensible heating was 
found to lead the stream function during maximum growth of 
the waves by a phase angle of 108-1149. This would put the 
Maximum heat input from the earth's surface to the 
atmosphere in a position upstream of the trough position and 
the zone of maximum cooling upstream of the ridge during 
amplification. 

In those instances when friction was not included in 
the modelling the vertical velocity wave led the stream 
function by approximately 102° indicating that maximun 
downward velocities in the mid-troposphere occurred between 
the ridge and the downstream trough position. When friction 
was added this phase lead was reduced by 4 or 5 degrees. 

Figures 4.10 and a. 41 show the curves for a wavelength 
of L,=7000 km. Such a wavelength is representative of long 
waves where the degree of instability is decidedly less than 
that of intermediate wavelengths near 4000 km. In both 
cases shown the initial phase angles were 90°. Figure 4.10 
illustrates model 4 while Fig. 4.11 pertains to the full 


model. 
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Fig. 4.10. Perturbation amplitudes and phase angles as functions of 
time for model 4 with lL = 7000 km, U = 20 m/s, UL = 10 m/s, Stee 90 
and « = 90°. Amplitudes in (a) are v (101m) é (°c) é CC) 

oO : A if ee 2 A : 


w,(10 * kPa/s) and h,(10 2 J/kg s)). Phase angles in (b) are in 
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degrees. 
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Fig. 4.11. Perturbation amplitudes and phase angles as functions of 


time for the full model with = = 90, Ly = 7000 km, U = 20 m/s, 


Oo 


Uy =O m/s, 6. = 90° and oe oe 90°. Amplitudes in (a) are (101m) , 


ts) 
6,(°C) , oncat ye w, (10-4 kPa/s) and He ClCie J/(kg s)). Phase angles 


in (b) are in degrees. 
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Several noteworthy differences are seen in these curves 
when they are compared with the previous resuits for L,=4000 
km. AS expected the perturbations now underwent a far 
lesser rate of growth than seen previously and tity moved at 
slower speeds as depicted by the shallower slope to the 
phase angle curves. As well, the rate of amplification of 
the waves now generally decreased with time and, in some 
instances, as illustrated in Fig. 4.10b, the waves 
experienced retrogressive tendencies, 

The first perturbation quantity to show a decrease in 
amplitude in Fig. 4.10a was the potential temperature. This 
occurred, according to the phase relations of Fig. 4.10b, 
when the potential temperature wave moved to within 20-25° 
of the stream function wave position. Then, as the 
potential temperature proceeded to overtake the stream 
function, the value of , also showed a reduction in 
amplitude. This is consistent with the concept of the 
thermal wave overtaking the height wave and thereby 
terminating amplification. In the latter part of the 72- 
hour period the potential temperature wave continued 
steadily moving away from the stream function which in turn 
began to retrogress slowly. 

Previously the perturbation field with L,=7000 km was 
observed to have greater instability and faster propagation 
rates once a meridional variation was inserted into the 
modelling. These features are now reflected in Figs. 4.11la 


and 4.11b where by t=72 hours the perturbation quantities 
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had all undergone greater amplification and had progressed 
further east than in the preceding case. The increased 
amplification is the result of the fact that the potential 
temperature never quite cvertook the strean anes wave 
and maintained a minimum lag of 22°, 

Changing the initial phase angles when L,=7000 km 
invoked much the same behavior in the phase angle and 
amplitude curves as was observed with the shorter wavelength 
of 4000 km. Greater amplification was found to occur when 


€,was taken to be positive (lag) and & negative (lead) 


6 
rather than reversing this arrangement. Of those numerical 
computations performed where sensible heating was included, 
the longer 7000 km waves showed the greatest amplification 
rates when €, was 180°. This arrangement implies that the 
thermal wave is initially exactly out of phase with the 
stream function wave and thus a warm trough and cold ridge 
structure is present. 

In most cases 72 hours was not enough time to establish 
any definite phase relationship between the static stability 
and stream function for the longer 7000 km wave. During 

A 
times of steady increases in ¥, the vertical velocity wave 
was found typically to precede the stream function wave by 
125° to 155°. This lead was shortened to as little as 60° 
when the stream function experienced weaker rates of growth 
or underwent slight damping. It should also be noted that 
the magnitude of the vertical velocity field was almost an 


order of magnitude less than that found at the shorter 
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wavelength L,=4000 km. The surface heating wave tended to 
remained downstream of the ridge position in the 
calculations with a leading phase angle of 120° to 130° with 
respect to the stream function wave. | 

At wavelengths longer than 7000 km the perturbation 
field continued its trend to slower growth and greater 
retrogressive motion. The 72-hour period of calculations 
was insufficient to determine many consistent phase 
relationships between the quantities although the behavior 


of the waves is somewhat similar to the L,=7000 km results. 


4.6 Blocking Ridge Development Over the North Pacific 
Blocking ridge activity, as interpreted by Rex (1950a), 
exhibits several distinguishing characteristics in its mid- 
tropospheric flow. These characteristics were defined as: 
(i) a Sharp transition from more or less zonal 
westerlies upstream to a meridional type of 
flow downstream must occur with a block; 
(ii) the basic westerly flow should have split 
into two branch streams with each branch 
transporting appreciable mass; and 
(iii) the two-jet pattern must exist with an 
extension of 459 of longitude or more 
downstream from the split and must persist for 
a minimal ten days. 
As an example of such ridging occurrences a detailed case 
study of an Atlantic block was presented in the Rex (1950a) 
investigation. This particular ridging lasted for more than 


20 days of which the first week could be described as being 


the growth time required to develop a recognizable block 
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from a basic zonal flow. 
In subsequent investigations Rex (1950b,1951) found 


that blocking showed a seasonal and spatial preference for 


its formation; preferring the winter months over that of the 


Summer and an ocean envircnment over that of a continental 
location. He also determined that blocking normally 
persisted for 2 to 3 weeks after its formation with very 
little overall motion during its entire lifetime. 

Whereas Rex employed daily synoptic charts in his 
studies, White and Clark (1975) used monthly mean data to 
investigate the occurrence of blocking activity over the 
North Pacific. This precluded the authors’ ability to 
determine specific growth or duration times for blocking and 
it could only he established that the blocking had occurred 
long enough to have affected the mean monthly pressure 
pattern in the mid-troposphere. In defining what was 
_considered as an oceanic block this study excluded all 
ridging waves which formed and were centered over the 
western part of North America as being a form of amplified 
mountain waves. 

Some of White and Clark's results have been mentioned 
in Chapter 1. AS shown in Fig. 4.12 the data revealed a 
preference for winter blocks which had a characteristic 
wavelength of 6000 to 7000 km centered over a longitudinal 
position of 160W to 170W. The blocking activity was also 
found rarely to exceed two months in length at any given 


period of the year and was apatially correlated with 
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Fig. 4.12. Seasonal and spatial variability of North Pacific blocking 


ridges as given by White and Clark (1975). The histograms were based 
on monthly mean data over a period from 1950 to 1970. The seasons 
were defined as being; winter (January-March), spring (April-June) , 


summer (July-September) and autumn (October-December) . 


patterns of ancmalous surface sensible heat flux (cooling 
under the ridge positicn). It was the authors! conclusion 
that Haltiner's (1967) predictions of slow moving or 
stationary unstable longwave growth explained some of these 
observations and, therefore, supplied theoretical support 
for their findings. However, upon closer look at the 
argument presented there would seem to have been some 
misinterpretation of Haltiner's results on their part. 


During their discussion White and Clark presented Fig 


82 


Tt 


eae pte ete 


Bae ou dog 46 oie sana “ he 
j { a Vr ’ 


i abit 


Gniaeotd aitige® : dsvon 3 veatdateas Reis bos ¥ 
bessd Srey vemeryoreta aAT hese Arad bain 04 v 
’ apioeade wel _ Be oze mort a 

, bamute ga) yalige , (savy) pari le a 


- (radnsest=v4803 209 oo 


seoeorr Tt a 


re ) pil “i oie 
“etter _s tued sldtenoe ebinnien avolsuoas 20 agyedtaq — 


Hy. 
ceceieaninel (21bAso8 edt new 9%. .inodketeng spbis edt rebay ‘Kas 


20 Yaivom wole Le enoisoiborg (Walty e* senldisd ond 
ves? ae save baatelgns dsuorng vvewpass sldszaay visnotsage 


i Bailaque -orehereds bas alla 


PAs, 99. 4008 xeobso 09H vteveron: sepabboss zkodd) 208 
"aden. gaed aved @? seek fisgoy axeds Sernees7q sesnupas. 
or deeds Go eeluros etyenitisn to sotasserqzetabake Pe 
Pr tasinnss saa fas pride aoleaugerh aied> exkanea ree 


4.13 as Haltiner's unstable phase velocity solutions for 


PHASE SPEED ( m sec”! ) 


ree ota = See Go bia To oa tO 
WAVE LENGTH (10° km_ ) 


Fig. 4.13.. The real (C,) and complex (C,) parts of the phase speed 
for baroclinic unstable waves as appeared in White and Clark (1975). 
The diabatic and adiabatic stability curves are denoted by Cr and 
Ci» respectively. Values are based on Haltiner's (1967) numerical 
results. The mean vertical wind shear and mean zonal wind were 


f , : A Ape 
given as uU, = 10 m/s and U = 30 m/s in White and Clark's discussion, 
however, Haltiner's theory utilized a mean zonal wind of only 20 m/s 


in obtaining these curves. 


U=30 m/s, U0 m/s and conditions with sensible heat 
transfer and without. In reality these curves are obtained 
for a mean zonal wind of only 20 m/s. It is also left to 
interpretation in the White and Clark study that the real 
portion of the phase velocity (C,) of Fig. 4.13 applies to 


both the adiabatic and diabatic cases. This is misleading 
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since the C, curve shown is valid only for the adiabatic 
case. The diabatic version of this curve for U=20 m/s gives 
a Stationary wave of approximately 7600 km which is some 300 
to 400 km shorter than the 8000 km found for the adiabatic 
curve. If U is taken to be 30 m/s as is used throughout the 
White and Clark discussion, the diabatic model as indicated 
in Table 4 gives a stationary unstable wave of L, =9000 km. 
This would then be well beyond the 8000 km value shown in 
Fig. 4.13. 

The behavior of Haltiner's amplified waves has been 
likened to that of Rossby-type nee. The Rossby dispersion 


relation of; 
1/2 
Zi = 27 aes (4.6.1) 


might therefore be used to predict the wavelength, L of a 


is 
stationary Rossby-type wave. Applying (4.6.1) with U=30 m/s 
and 6 (45N)=1.61X10~-1! m~1s-1 gives a predicted L, value of 
8577 km. If U=20 n/s is substituted into the relation the 
stationary wavelength diminishes to 7003 km. Neither of 
these values agrees with a value of 7600 kn claimed by White 
and Clark using fabeenty mand U=30 m/s. The apparent 
confusion surrounding White and Clark's argument makes it 
necessary to reconsider the theoretical predictions and 
their applicability to the characteristics of North Pacific 
blocking development. 


In order to obtain values for the mean winter zonal 


wind and vertical wind shear to be employed in the 
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Fig. 4.14. The zonal average of the mean westerly winds for January 
at 80 and 40 kPa over the Northern Hemisphere as given by Gates (1975). 
The open circles are zonal averages of the observed mean wind (for 
December, January, February) based on the data of Crutcher (1961). 

The solid curves represent the corresponding zonal averages of the 


observed geostrophic winds as based on the data of Crutcher and 


Meserve (1970). 


computations, use was made of the curves of Fig. 4.14 which 
are taken from Gates (1975). These curves are based an 
available observed wind data in the Northern Hemisphere and 
those mean geostrophic zonal winds which may be calculated 
from observed geopotential fields. As shown, the zonal 
average of the 80 kPa westerly wind had a definite maxinun 
value of 8 m/s at approximately 4ON to 45N. The observed 


westerly wind at 40 kPa reached its maximum of just over ZA 
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m/S at 30N to 35N although the geostrophic wind was somewhat 
stronger at 28 m/s. In terms of vertical wind shear across 
the two isobaric levels, the greatest difference in the 
observed mean westerlies was found to occur near 35N and was 
about 16 m/s. For the geostrophic zonal wind the maximum 
difference was slightly higher being 20 m/s at about the 


same latitude, 


aR tote aft: oie Ay ew) 
U Zonal wind (m/s) 
Fig. 4.15. Vertical profile of the observed mean zonal wind for the 
period December to February at 45N as given by Déds (1962). Values 
are based on the data of Crutcher (1959). Pressure (ordinate) is 


given as a non-dimensional variable P = p/P. where Pee 90 kPa. 


Figure 4.15 taken from DOd6s (1962) is also based on 
averaged wind data. It shows the observed zonal wind 


profile for the winter period December to February at a 
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latitude of 45N. The profile displays an almost constant 
wind shear through the mid+layers of the troposphere which 
translates to roughly 12 m/s over the 36 kPa interval 
between P=0.8 (72 kPa) and P=0.4 (36 kPa). 

Based on Figs. 4.14 and 4.15 a representative mean 60 
kPa winter value of U at 45N would be approximately 13 or 14 
m/s. Since observed year-to-year changes in the 
autumn/winter mean wind in the lower troposphere are 
typically in the order of 25% (White and Clark,1975) an 
assumed range of values for U was taken to be 14 = (.25)X(14) 
M/S, giving a variation from 10.5 to 17.5 m/s. The 
corresponding mean vertical wind shear in the nid- 
troposphere was taken as being about 7 m/s over a 20 kPa 
layer. Normally the mean wind shear in the atmosphere at 
mid-latitudes does not experience large fluctuations. White 
and Clark have established that the year-to-year variability 
in the winter value of the vertical shear for the 
geostrophic wind between 100 and 70 kPa over the North 
Pacific is only 5%. Since it is not inconceivable that this 
also applies to the mid-troposphere a range of values was 
not assigned to U_. | 

It is worth noting that by using U_=7 m/s the implied 
mean westerly wind at 50 kPa is close to 17 m/s. If a 25% 
variation is again instituted the upper estimate becomes 21 
m/s which falls significantly short of White and Clark's 
estimate for the mean winter westerlies of 30 m/s at this 


level. General circulation studies for January, such as 
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those of Gates(1975) and others, typically show a simulated 
local mean westerly wind maximum off the east coast of Asia 
at mid-latitudes which may extend well into the North 
Pacific. If correct, this implies that the winter value of 
U may be slightly higher in this region than that derived 
from the zonally averaged wind data. It is not likely, 
however, that the increase in the mean wind over the entire 
North Pacific is enough to account for the rather large 
difference between the averaged values and that of 30 m/s. 
It would appear from previous studies by other authors that 
the lower value of near 20 m/s is more commonly used when an 
estimate for the mid-latitude 50 kPa mean wind is required. 

Earlier it was seen that the modelling is not 
critically responsive at longer wavelengths to changes in 
the mean static stability. Therefore, o was held constant 
at 12°C in the computations. Discussion presented in White 
and Clark's study indicated that the surface sensible heat 
exchange does not vary significantly during the winter 
months. Hence, k, was also given a constant value of 10-2 
J/(kg s °C) during the calculations. 

The adiabatic version (model 1) did not predict the 
possibility of stationary unstable wave growth when U_=7 m/s 
and U=10.5 to 17.5 m/s were used in the calculations. This 
was not the case when sensible heating was incorporated. 
Table 8 summarizes the phase velocity solutions of model 2 
for those wavelengths with which the possibility of a 


stationary unstable existed. From this table it is 
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Table 8. Unstable wavelengths (L.) and the corresponding e-folding 
times (t.) for model 2 (diabatic, frictionless flow) using U, = 7 m/s 
and ke = 10 2 J/ (kg S Be) over a relative phase velocity range of 


approximately -10.5 to -17.5 m/s. 


L. (km) t, (days) Cp, 7 U (m/s) 
82 eS nhl Say Wee eine See 3 ee 
5500 6.5 ~10.1 
5600 6.9 -10.6 
5700 7.3 “11.1 
5800 7.8 “11.6 
5900 8.2 -12.1 
6000 8.8 “12.7 
6100 9.3 -13.2 
6200 9.7 -13.7 
6300 10.4 -14.3 
6400 10.9 -14.8 
6500 WS -15.4 
6600 12.0 -16.0 
6700 12.7 -16.6 
6800 13.3 “17.1 
6900 14.0 “17.7 
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conceivable that a wave within a wavelength band of 5500 to 
6800 km could undergo stationary amplification. The upper 
end of this waveband would require a slightly Alone auc 
westerly wind and longer growth times while aries 
wavelengths would give fairly rapid development with a 
weaker zonal flow. In general the waveband was centered 
about a value of L,=6200 km which had a relative wave speed 
of about 14 m/s and an e-folding time of approximately 10 
days. These values are in good agreement with 
characteristic growth times of blocking (1 to 2 weeks), the 
previous estimate of the mean westerlies and observations of 
winter blocking ridge activity over the North Pacific, 
Model 2, however, includes only sensibie heating effects. 
The question therefore arises as to whether incorporation of 
friction in the modelling might also give similar 
predictions. 

Table 9 reviews the results of model 3 for U,=7 m/s. 
The addition of friction to the adiabatic flow produced only 
a very narrow band of possibilities for stationary wave 
growth, ranging from 6000 to 6300 km. A look at the e- 
folding times shows that, except for L,=6000 km, the growth 
rates of these waves were very slow. It is, therefore, 
unlikely that friction alone could result in stationary wave 
amplification within the model of the same time and space 
scales as atmospheric blocking. 

Calculations were dea done for model 4 which 


incorporates sensible heat exchange and friction 
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Table 9. Unstable wavelengths (L.) and the corresponding e- 
folding times (t.) for model 3 (adiabatic flow with friction) 


using U. = 7 m/s and k. = 4 x 10 & s ! over a relative phase 


T 
velocity range of approximately -10.5 to -17.5 m/s. 


L,. (km) t, (days) Caz ~ U (m/s) 
5900 10.0 -10.1 
6000 14.2 = OeL 
6100 22,5 “11.3 
6200 42.3 -11.9 
6300 166 -12.6 
6400 > SUSior 


simultaneously. The waveband of interest was determined as 
being from 5500 to 6200 km (Table 10). In almost: all cases 
within this waveband the mean westerly wind was necessarily 
slightly less than its long term mean of 14 m/s for 
stationary wave Aeeonnanes The growth rates of these 
waves, as indicated in Table 10, were a little too slow to 
describe blocking development and the central value of the 
range of wavelengths was a rather short 5800 km, Such a 
wave would require a mean westerly wind of only 12 n/s to 
become quasistationary, but almost a month to increase its 


amplitude by a factor of e. The addition Of frictional 
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Table 10. Unstable wavelengths (LY) and the corresponding e-folding 


times (t.) for model 4 using U. = 7 m/s. Surface heating and 


friction constants used were k. = 10 2 J/(kg m SC) and k =4x10 & 


so. respectively. The relative phase velocity range is approximately 


-10.5 to -17.5 m/s. 


L. (km) Rs (days) CR3 7 U (m/s) 
5400 14.4 =1023 
5500 16.6 -10.8 
5600 1950 yb 
5700 23.9 -11.8 
5800 73} | =e 
5900 38.8 =1257 
6000 55 —3.3 
6100 86.4 =ae0 
6200 228 -14.3 
6300 - -14.8 


effects to the model has to a certain extent offset some of 
the congruity between the predictions of model 2 and 
observed characteristics of North Pacific blocking. 

It should be cautioned that the value chosen for U, in 
these calculations was critical. For instance, Table i1 is 
taken from model 4 under the same conditions as before 


except that now U.=8 m/s. The range of possible stationary 
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Table 11. Unstable wavelengths (L.) and the corresponding e-folding 
times (t.) for model 4 using k, = 10 2 J/(kg m/s) and k. =h x 10 & 
s 1, Vertical wind shear has been increased to UL = 8 m/s in the 


calculations. The relative phase velocity range remains unchanged 


atsrougnly <10.5 to -17.5 m/s. 


L. (km) ce. (days) CR3 U (m/s) 
5500 6.2 -10.4 
5600 8.7 -10.9 
5700 9.7 -11.3 
5800 10.8 -11.8 
5900 122 -12.3 
6000 14.0 -12.8 
6100 iy G8) 1 Bh) 
6200 1857 -13.9 
6300 2203) “14.4 
6400 27.4 -14.9 
6500 34.2 -15.5 
6600 45.7 -16.0 
6700 68.6 -16.6 
6800 114 -17.2 
6900 «42k -17.8 
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unstable waves was found to exist from 5600 to 6800 km. The 
central value of near 62C0O km had an e-folding time of 
roughly 19 days and required a mean zonal wind of 14 n/s to 
Maintain a quasistationary position. The ero Taunt growth 
time is slightly greater than found for atmospheric 
blocking, but ctherwise the values are in agreement once 
more with observations of winter blocking over the North 
Pacific and the estimate of the mean westerly wind speed. 
Therefore, if friction is included in the modelling in 
addition to sensible heating, a 1 m/s increase in the value 
of U, was required to overcome its influence. This might 
suggest that atmospheric blocking occurs in the presence of 
friction if the mean vertical wind shear in the nid- 
troposphere is slightly above its long term mean value. 
Without the presence of sensible heating, however, 
stationary wave growth as predicted by the modelling is not 
likely to occur. This may explain the observed decrease in 
atmospheric blocking activity over the North Pacific during 
the summer months since at this time sensible heat exchange 
between the ocean and atmosphere is greatly reduced. 

In studies where the mid-latitude standing wave problen 
has been investigated by means of linear steady-state 
models, it has been found that certain meridional 
wavelengths for the harmonic perturbations give more 
reasonable results than others when substituted into the 
equations. Such was the case in the present study for 


computations based on the full model version of the 
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equations. Table 12 summarizes the results for three 
different meridional wavelengths: Ly =60°, 90° and 120°. The 
value of U, was again taken to be 7 m/s. With ae the 
full model does not predict the possibility of a stationary 
unstable eae within the 10.5 to 17.5 m/s range selected for 
the winter value of U. When the meridional wavelength is 
changed to 90° a small band of unstable stationary waves was 
possible between L,=6400 and 7000 km. However, the growth 
rates for the waves within this waveband were very slow, 
being of the order of one to several months. The last 
value, bya 200s yielded possible stationary wave 
amplification between approximately 5900 and 6700 km. Once 
again the e-folding tires were fairly large, although at the 
shorter end of the waveband t, had decreased to about 3 
weeks. 

In general the full model, when used with U;=7 m/s and 
L, =6 0° or 90°, did not satisfactorily describe winter 
blocking activity over the North Pacific. With L,=120° the 
full model predictions were similar to those of model 4 
except in one significant regard, this being the value of 
the central wavelength within the waveband of possible 
stationary unstable amplification. In the last column of 
Table 12 this central wavelength is in the vicinity of 6300 
km which is some 400 to 500 km longer than the central 
wavelength predicted by model 4 (Table 10), yet the e- 
folding times are comparable . The longer 6300 km value 


lies in the middle of the observed mean wavelength range of 
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Table 12. Unstable wave solutions for the full model using 


U, = 7 m/s, k, = 102 J/(kg s °C) and k, = 4x 106s, The y-axis 


t 


wavelengths (L.) are given in degrees of latitude with the origin 
y = 0 at 45N. The x-axis wavelength (L.) interval of interest is 


between 5500 and 7500 km. Relative phase velocities (c - U) are 


R3 


given in units of m/s with e-folding times (t.) expressed in units 


of days. 
Les 60° Las 90° 

L, (km) ee Croaanl mn easel 

5500 8.4 -6.1 VOn3 -8.0 .0 
5600 | 9.0 -6.2 13 -8.3 A 
5700 9.6 ~6.4 25 -8.6 a 
5800 10.3 -6.5 13.9 -8.8 al 
5900 11.0 -6.7 15.3 -9.1 A 
6000 11.9 -6.8 |W hes) -9.4 8 
6100 12.8 -6.9 19.4 -9.6 ri 
6200 ot -7.1 22.0 -9.9 as 
6300 14.9 -7.3 2552 =30.2 .8 
6400 16.1 -7.5 29.5 -10.5 2 
6500 17.6 -7.7 35.2 -10.7 .6 
6600 19.6 -7.8 43,4 -11.0 9 
6700 2\a7 -7.9 56.0 -11.3 aS: 
6800 23.6 -8.1 78.3 --11.6 G 
6900 26.5 -8.2 114 -11.8 

7000 30.0 -8.4 258 -12.1 

7100 34.4 -8.5 - -12.4 

7200 hO.2 -8.6 

7300 46.4 -8.7 

7400 56.8 -8.9 

7500 hed) -9.0 
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Table 13. Unstable wave solutions for. the full model OLth the 

mean vertical wind shear, UL = 8 m/s. The respective sapees 
constants for sensible heating and friction are k, = 10 2 J/(kg s ec} 
and kK. = 4x 106s 1, The y-axis dependence is taken as being 120 
degrees of latitude. Values shown are for a relative velocity 


range of approximately -10.5 to -17.5 m/s. 


Ly (km) t. (days) CR3 7 U (m/s) 
6000 9.6 -10.4 
6100 10.6 -10.7 
6200 11.6 11.) 
6300 12.9 “11.4 
6400 14.2 -11.8 
6500 16.0 -12.1 
6600 hes -12.5 
6700 20.6 =12.9 
6800 23.6 =13%2 
6900 27.6 -13.6. 
7000 Bake -14.0 
7100 39.6 -14,3 
7200 | 51.0 -14.7 
7300 6722 15.1 
7400 97.4 =o 
7500 173 -15.8 
7600 700 -16.2 
7700 - -16.6 
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North Pacific blocking which is 6000 to 7000 km according to 


Fig. 4.12. 

Recomputed predictions for the full model with L,=1209 
and U.=8 m/s showed an improvement in the growth rates as 
was seen with model 4. AS given in Table 13, all waves in 
the 6000 to 7000 km range lie within the 10.5 to 17.5 m/s 
range for U and could conceivably undergo stationary 


amplification with e-folding times ranging from 10 days at 


the lower end to roughly one month at the upper extreme. 
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CHAPTER 5 
SUMMARY AND CONCLUSIONS 


In this investigation a linear two-level baroclinic 
model was formulated which incorporated a self-determining 
variable static stability. The attraction of this model was 
that the adiabatic-frictionless form of the equations 
contained certain invariant energy and vorticity properties 
which were desirable. Sensible heat exchange and frictional 
effects were then introduced in the equations as 
perturbation features within the steady-state current flow. 

Harmonic solutions, independent of latitude or assigned 
a simple meridional wavelength variation, were assumed for 
the perturbation quantities. A cubic frequency equation was 
obtained and sclved for varying conditions of current flow. 
Analysis of the resulting phase velocity response curves 
showed that: 

(i) adiabatic-frictionless flow yielded a band of 

unstable waves with maximum instabil wey sat 
shorter wavelengths; such flow also exhibited 
a longwave cutoff, above which wave 
amplification was not possible. 

(ii) sensible heat exchange reduced the 

instability of the unstable waveband at short 


and intermediate wavelengths while increasing 
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the instability of longer waves. 

(iii) friction diminished the instability of the 
perturbation field at all wavelengths except 
in two narrow regions found at the extremes of 
the unstable waveband. 

(iv) simultaneous inclusion of friction and 
sensible heat exchange created the greatest 
reduction in the amount of instability 
available for short to medium length waves, 
but resulted in a longwave cutoff of longer 
proportions than that found with adiabatic- 
frictionless flow; coincident with this the 
propagation speeds of the longer unstable 
waves were significantly slower. 

(v) the instability and wave speeds of the long 
waves increased as the meridional wavelength 
of the perturbations decreased. 

The stability of the perturbation field was also 
subject to changes in the steady-state flow. Most notable 
were the effects of using larger vertical wind shear values 
and shifting the x-axis to higher latitudes both of which 
created increased instability in tHe waves. A reduction in 
the atmospheric mean stability produced more instability at 
shorter wavelengths, but had appreciably less effect at 
longer wavelengths. , 

An initial-value technique was employed to determine 


the wave amplitude and phase characteristics of the various 
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perturbation quantities. Numerical solutions supported the 
findings of the stability analysis and indicated lesser wave 
amplification and slower wave speeds as wavelengths 
increased. It was also established that the thermal wave 
necessarily lags the stream function during unstable wave 
growth and that, conversely, the static stability leads the 
stream function wave. Maximum downward vertical velocity 
and the location of maximum surface heat exchange to the 
atmosphere were found between the ridge and the downstream 
trough position. 

It had been suggested in a previous study by White and 
Clark (1975) that regional North Pacific blocking activity 
owed its existence to its marine environment and the 
presence of sensible heat flux. The present study 
investigated the space and time scales of the unstable 
baroclinic long waves which were predicted by the theory. 
Using representative mean winter values for the zonal wind 
and vertical wind shear, the numerical values indicated that 
stationary unstable wave amplification having similar 
Characteristic growth rates and wavelengths as North Pacific 
blocking ridges was possible when sensible heating was 
present. Agreement between the numerical results and 
observed values was improved if: 

(a) the vertical wind shear in the mid 

troposphere was marginally above its long tern 
mean value; 


(b) the perturbation field was given a single 
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meridional wavelength of 120° of latitude. 
In the absence of sensible heat exchange stationary wave 
growth was possible, but not realistic because of the 
extremely slow growth rates involved. Certainly such 
amplificatjon would not occur on the same space or time 
scales as observed with blocking ridges. 

The position of naxinun heat exchange within the 
modelled atmosphere supported the spatial correlations 
between blocking and anomalous heat flux distributions over 
the North Pacific ocean found by White and Clark (1975). 
They found that anomalous surface sensible heating was 
negatively correlated with the position of a blocking ridge 
formation. The initial-value technique indicated that at a 
wavelength of 7000 km maximum perturbation heating in the 
surface layer of the atmosphere occurred some 50° upstream 
of the trough position during development, thereby putting 
maximum cooling just downstream of the ridge. 

The concept of blocking ridge activity as a form of 
baroclinic instability would not seem unreasonable from the 
results which were found in this study. However, the theory 
presented can only cast inferences on the development of 
blocking formations. Linear perturbation theory can say 
very little of the ridge's ability for sustained existence 
without apparent further amplification. The model used in 
this study excluded a major mechanism for heating in the 
atmosphere — that of latent heat release. Further 


investigation is required to determine in what way this heat 
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source would influence the numerical predictions. 
Orographic and land/sea constraints were aiso neglected in 
the modelling. Kikuchi (1971) found in his numerical 
studies that mountains and land-sea distribution played a 
major role in the longitudinal position of the blocking 
action, but were not necessary for the creation of blocking 
waveforms. In addition, his study indicated that 
nhumerically-simulated blccking ridges had a tendency for 
longer duration if terrain effects were included in the 
modelling. Additional theoretical approaches are required, 
however, to substantiate the role that orographic effects 


play in the scheme of blocking formations. 
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